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AFPEHDIX A 

I'ELLCEIT’S ADVANTAGE (MULTIPLEX ADVANTAGE) 


A.i TUEOUy 

Fellgett's Advantage is the name of a factor which shows the improvement 
obtained when the spoctrum oF thermal radiation is moasurod by a Fourier traiis- 
l.'urm spuctrumuter^ ratlier clian by a scanning filer spectrometer. 

This difference in performance only occurs when there is no amplifier 
ahead of the spectrometer. Then Fellgett'K Advantage i.s a measure of spec- 
trometer performance when all of the input power is applied to a single 
detector in the autocorrelator compared to the performance obtained when a 
narrow, bandpass filter scans the spectnaa. Such a filter wastes the spec- 
trum power outside the filter passband. 

Fellgett's Advantage is found as follows. Let 


D ■ total system bandwidth 

Tsys “ system noise temperature (A-1) 

N - k TjjYgB - noise power (A-2) 

T ■ total observation tine 

e 

M "■ number of frequency bins resolved 

Then in the case of the scanning filter 

b ■ B/M ■ radiometer bandwidth using a scanning filter 

T " T/M ■ time available to measure temperature of each frequency bin 
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(A-3) 


Tn ulie correlator case 
b M B 
T 

T * u " available for measurement of each correlation sample* 

Since is complex, a real end an imaginary sample are 
tukuo. Tim specurum generated by Fourier transformation has 
the same number of samples. However, the samples are real. 
Therefore, the variance of spectnun bln is the sum of the 
variances of the ronl and Imaginary samplos, 


^'''coiin 
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(A-/.) 


Fellgett's Advantage, F, relates T 


FILTER ^COR 


AT 


CORK 


— X 

M ^FILTER 


(A’"5) 


Equation (A-5) states Chat if we cannot increase the power at our dis- 
posal by amplification, Che noise level of a scanning filter spectrometer is 
greater than that of a Fourier spectrometer. 

However, if amplification is available, power is available to feed a 
channelizing filter bank or a delay line with a set of taps. These devices 
will iJive a temperature resolution which does not depend on H. 

* With no restrlcLluiv on the number of taps, we can measure i|>(t) for both 
positive and negative values of t. Now is symmetric and ♦ (t) imaginary 

is antisymmetric. So the corresponding positive and negative real correlation^ 
sattpleS, can be added together. The imaginary samples paired with the 

real samples can be added after reversal of the sign of one sample. In both 
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cases Che noise adds noncoherenCly and che sigtuil coherently. This reduces 
the temperature resolution by 2"l/2. This compensates for the 2 Increase in 
the Fourier spectrometer case caused by the fact that is complex, but 
the power-density spectrum is real. 

The T * 0 tap is a special case. The <^(0) jjiiaginary s«“plc must be zero. 
So It is not used in computing the transform. This produces a two-to-one 
noise reduction consistent with the noise reduction obtained with the pairs 
of caps wljcrc the dulnys nre (t 0) . 

A. 2 PERFORMANCE CALCULATIONS 

• As examples of the relative performance of these devices let us consider 
three correlator/spectrometers operating between the water vapor line and the 
©2 line complex, say 23-53 GHz. 

One device will be a crystal video system using a Hughes Electron 
Dynamics GaAs diode,' a second system uses an evacuated barretter, a third sys- 
tem uses a TRG mixer and a Narda GaAs FET IF amplifier. The amplifier has an 
8-18 GHz passband. 

This amplifier has a noise figure of about 9 db according to the NARDA 
ads.'^^. The Trk mixer appears to have a single sideband noise figure of about 
6 Tills supcriict must upurzitu in u single sidebun^d mode. Thu double 

sideband mode would put a 16 GHz hole in the middle of the reception band. 

The total noise figure is then about 15 db. 

A TSS of 6 X 10*" watts will be used for the evacuated barretter . 

_9 

Figure A-T shows a TSS oF -'>5 dbm or x 10 wattii Is a good 

uuuibuc I'ur Llie dludu purl'urtuuncu. 

Performance calculations for the three cases follow.- A radiometric per- 
formance parameter of K > 1, and an integration time of one second will be 
used in all three cases. The result is a representative temperature resolu- 
tion which is useful for comparison purposes. It also will be assumed that 
1024 frequency bins are desired. It then tehee 17.07 minutes to measure the 
spectrum. 


Microwaves, p. 7, April 1980. 

Gardlasmenos , "Planar Devices Make Production Practical," Microwave 
Systems News , April 1979, pp. 15-21. 

^^^Long and Butterworth, "New Technique for Microwave Radlometry," IEEE 
TRANS MTT, Sept 1963, pp. 389-396. 
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FREQ, OHZ 

NOTE! oiAS cunnENT in a I'AnAwnrcn. tss is «i» to g mi iiuttch 
WHEN TUNUU MOH NAHHOW UAND RESPONSE « IOGH 2 ). 


FIGURE A-1. TSS VS FREQUENCY 

SUl’liUUJiT 

W - 16 GHz 
NT - 15 db 

^SYS " * 290“K - 9170. 61 *K 


FILTER BANK OR TAPPED 

DELAY LINE AT’ - 

- 7=1 - 0.070*K 



/i; 

SCANNING DELAY 

AT - 

✓5o48 X AT* - 3.28®K 

SCANNING FILTER 

AT - 

1024 X AT’ - 71.68-K 


BARRETTER 


T-vc " -i ^ T7? - 218.98X AT’ - 0.00173*K 

X.37 X X 2 X 10^ • 


FILTER BANK OR TAPPED LINE NOT POSSIBLE 
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SCANNING DELAY 

* 


AT ■ /2048 X 1.73 x lO”"* - 0.08*K 

SCANNING FILTER 


AT - 1024 X 1.73 X 10"^ - 1.77*K 

CRYSTAL VIDEO 


* 

3.16 X 10’*'* W 


- 11,533*K, AT’ - 0.0815*K 

1.37 X 10’^^ X 2 X 

10^° 

FILTER BANK OR TAPPED DELAY 

LINE 

NOT POSSIBLE 

SCANNING DELAY 


AT - /2W X 8.15 X 10"^ - 3.69’K 

S0ANNJ.N(i FILTHH 


AT « 1024 X 8.15 X lO"^ - 83.46*K 


It may be seen that the acnnnlriQ filter mechanization is a very poor per- 
furiuur. i{oy(.*ver, cite Kourleir tcunsl'onn .syscem with a scanning delay is a 
remsonable system mechanization, if sinipllcity and cost are important. In 
cases where the ultimate in performance is not needed, it is a viable candidate. 
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APPENDIX B 

GENEliALXZEU XNTEKFEKDHETER NAPPING BEHAVIOR 


The basic properties of correlating iuterferorH^ters as mappers should be 
addressed. The fact that the system is linear in its response to power or 
intensity will be used in several key derivations. To show under what condi- 
tions such linear rusponso is obtained » a calculation of the mapping response 
of a correlation radiometer in the general case follows. 

The interferometer system is indicated below. 

Let F(x, t) be the electric field distribution across the source. 

Assume it to be a stationary random process. Truncate this process to form 
F(T, Xf t) so that a function with a Fourier transform is produced. This 
fiuid is complex and analytic. The sccund property arises because the real 
and Imaginary components of ? are assumed to be Hilbert transforms. This 
guarantees that F is a causal, oignal. 

Thu Fourier transform of V(x, t) is 


F(T^, X, V) 


. fnr. 


X. t) e 


2irivt 


(B-l) 
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FIGURE B-1. CORRELAtING INTERFEROMETER 
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wiluco V lti:u(|uunuy, u « tXino, X > uuurcii coordinabe and hae tynbol Indlcabas 
cransform funccion. 

Zn order bo charmebarize eh« antannaB, lab G]^<x) and C 2 (x) ba tha aaipll- 
buda diffracbion pabborna oF, A] nnd A 2 > Bobh pabbornR nru in a polar coordi- 
imbw NyHbOiu wIlIi url^i.iii ub 0 In Pigucc B-1. i'hon blto oubpub of bha firab 
anbanna la 


Vj^CT, X, v) » v)>»'VCT, X, b) 


S lilt I larl y 


vj(T, X, v) » g|(x, v)*V*(T, X, b) (n-2) 


Tha producb of bhaae bwo signal produces bhe video signal output with 
blio form 


V^(T, X, \») vJ(T, X, V) 



V) C2^X2» 



V) V*(T, X2* 



(B-3) 


where ** indicates a two^diaensional convolution (* ■ one-dimensional 
convolution) . 

Divide both sides of Equation (B-3) by 2T. Since gJ and G« do nob depend 
on X, the stabisbics are stationary. In bha llmib ^ 

B(x, v) *2' (B-4) 

t-m. 


A ^ 

where P(x, v) is the complex spectral density of the output and TCxj^, X 2 » v) 
is the mutual spectral density of the radiation in the plane of the antennas. 
(The Van Cibber b-Zernike theorem describes how these coherence functions 
propagate from bhe object scene to the antenna plane.) By definition 




V(Xj^, \>) V(x 2 , v) 
2T 


11m 


(B-5) 


2T 
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Equation (B-4) is the gtnaral* aolution for tht intarfaronattr rosponaa 
CO an axCtndtdi, parclalXy cobac«ne» polychromatic aourca* 

for cha qumaiHionochroauiCic caaa, which holda for moat practical aya- 
Caaat thara ara two condiciona. Thaaa ara 


Av a: V 


Itl < ^ and i a? ^ (B-6) 

s 

That la, the spectral width, Av, must be small compared to the mean fre» 
quency, 7; also all system path differences must be small compared to c/Av> 

For this case the mutual coherence function becomes 

r(Xj^, Xjj, t) - r(xj^, X2» 0) (B-7) 


and the mutual power cross-spectral density is 

A A 


rCxj^, Xg, t) - r(xj_, Xj, 0) fi(v - ^ 


(B-8) 


Substituting (B-8) into (B-4) and evaluating at t 0 

P(x) - [G^(x) GgCx)] ** r(x^, X 2 , 0) (B-9) 


G;i^(x) and G 2 (x) have become antenna patterns of and A 2 at frequency v 
This is the suluticn fur the output response to a partially coherent, quasl*^ 
monochromatic, extended source distribution* 

Let us look at the limiting cases of coherent and noncoherent extended 
aourcas. For a coherent source 


0) - V(x^) 


(B-10) 


Beran and G. Parrent, Theory of Partial Coherence, Wiley & Sons, 1964, 


Sec. 4.3 and 4.4, pp. 54 and 5 
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Then (B>9) becomes 





I 

P(x) - [G 3 ^(x)*V(x^)] [G 2 (x 2 >*V*(x 2 )] 


(B-U) 


This equation shows that for mapping extended objects which are 
coherently illuminated (or radiate coherently across their surface) the Sys- 
tem is essentially nonlinear. Unscrambling such a nonlinear signal to, form 
an image appears to be a hopeless undertaking. 

In the l imit, wIhtc V(x) « ‘ctS(x - x.) for n point ftuuruu (B-11) reduces 
to 


P(x) - a.^U) a^(x) (B-12) 


Thus, for poi,iit sources, the system secs a power diffraction pattern 
which is the product of the patterns of the two antennas.' This appears 
intuitively to be correct. 

(•'or seir-l.umiiiouii Mourcos wliiuh are essuntiaiiy noncoherent 




The subscript "0" denotes that this is the distribution at the object 
source. 

Substituting (B-13) in (B-9) and integrating over the source, we get 


P(x) - y e P(j(x’) [G^ix - x') G*(x - x’)] dx' (B-14) 

•'source 


Where Pq(x’) is the distribution op power over tlics surfar.e uT the source 
Equation CB-14) describes a system which is linear in power and wlfich has a 
complex impulse response K(x) given by 

K(x) - C^(x) gJ(x 2 > (B-15) 


B-5 
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I'lion cfin be wcittett 


P<x) - JVq(x)*K(x) <B-16) 


Thus the system under consideretlon, although nonlinear for mapping 
extended coherent fjurces, is certainly linear for mapping extended incoherent 
sources such os the radio^^sky or rod iome trie-earth. 


Since and G 2 are defined relative to the origin 0, we can set the 
antenna phase centers at that point. Then 


! Gj^(x) - g’^(x) 

G2<x) « g '^( x ) 


(B-17) 


where g],(x) and G 2 (x) are the radiation patterns of A 3 , and A 2 relative to 
their own phase centers and ^ Is the separation be cwoon these phase centers. 
TiK*ruii>re, Lite iiiipuJiK' responiu* can be written In the I'uriit 


K - G^(x) + G^Cx) (B-18) 


hot Ax and A 2 itavu real» syituuetclc weightlug l'uuctiona» Ax(h) and A 2 (n)» 
where n ^he untcnnu aperture coordinate. Then and G 2 will also be real 
and symmetric. For this special case. > 


•^ssmuetrlc'*) ‘ ^Gi^) cos (B-19) 


In this special case the interferometer impulse response is simply the 
product of the amplitude patterns of the two antennas and the interference 
factor resulting from the phase center separation. 

By the cuavulutluh theureiu tlie Fourier transform of (B-16) is 


P(x) - P„<x) * K(x) 

if tj . 

S(F) - 0(F) • MTF(F) (B-20) 
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where 0,(F ) a nd S(F) are the complex spatial frequency spectrums of image and 
objectt’ MTp(F) Is the modulation transfer function of the intcrferomctor . 

Taking the Fourier transform of (B-15) 


K(x) - G (x) 

J I 

MTF(F) • A^(n) 



* A2(ri) 


(B-21) 


The second two transforms are the Fraunhoffer diffraction relation 
between aperture excitation and the resulting radiated field. 

Braceweil gets rid of the characteristic folding Inversion due to the 
convolution by using the mirror Image of K(x) . Gall It K(x). Then the bottom 
relation In (B-21) becomes an autocorrelation. Thus 


MTF(F) - Aj^(n) * A2(n) 


(B-22) 


where 

K(x) - Gj^(x) • G 2 <x) 


(B-23) 


It sliuuld be obvious that K(x) Is the antenna power pattern of the 
Interferometer. 
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APPENDIX C 

MICROWAVE RADIOMETER IMAGE STATISTICS 

I 


C.l SPECTRUM SELECTION USING RESULTS DUE TO MANDELBROT 

Mandelbrot's computer generated Brownian landscapes on pp. 210-215 of his 
book show how changing and Image parameter, H, changes the appearance of the 
terrain images. They show how the correct value of H Is needed In order to 
produce the subjective impression that one Is looking at an "image" and not 
noise. Unfortunately, these figures will not xerox. However, the computer 
genera Ced Brownian trails shown below could be reproduced. These trails show 
the change In long range coherence due to changing the parameter, H. 

Mandelbrot concludes that in order for his Brownian landscapes to have 
the proporlvioK of ronl laiulncaiiuH the redundancy pacamutur, U, Hiioiild have a 
value of 3/A. This corresponds to a spectrum of the form, 


G(£) « sf 


• 211-1 


» Sf 


•5/2 


The tcanKrorm-mate iiuLocorrclation expression Is given In Elgure C-2 along 
with a plot of this functlon.i 


C.2 CALCULATION OF REDUNDANCY PARAMETER FROM IMAGE C0MP1U2SSI0N DATA 

In order to calculate the properties of the various imaging systems, we 
need to use a power spectrum for the radiometer Images which has the right 
redundancy. 

Assume that all images (visual, IR, or microwave) have about the same 
statistical properties, because they originate from the same object scene. 
Then we can estimate li from knowledge of propertlas of images in cofiipressed 
format. Such compression is carried out in order to send Images with reduced 


Mandelbrot, Fractals — Form, Chance and Dimension, 'W. Freeman & Co., 
1977. " 

/ r?-, 

iv-' ' 
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The Flpire lo the left coiiithulef in 
iilmplc o( a iialifiically icif ilmilar 
fractal curve with the dlmennon 0* 
\ /0.9000* l.n I U lit cwxdtnaie («m* 
ilwM are independent fractional Urowni* 
M funciioni of eaponenl HaOtSOOO, 
which ImpUea that Mif intertccilont. 
thou|h not iirictly prohibitad. arc greatly 
dheoMraged by forctny the curve to be 
eery pentitent. Thmkmi of complicated 
ewrvea at the aupeiimpotltluni uf large, 
raar M ura. and tmall convolutiont, It may 
ha mM that in the cate of hiph pcniii> 
•aaa and diaternion claac to I, aaiall 


The Fi|urc to Ihc r{|ht utet the tame 
tencraiini program and the tame ptcudo 
random Htd but a dimenaion increaicd 
to 0 • I /0.7000* 1 .4388. Without hay. 
ing to change the ihape of the convolu. ^ 
iloni which add up to form thlt trail we 
have incrcaied the relailra importande of 
Ihc tmall onci, and to a Ittser eitcni, of 
the medium onei, Previoutiy Invlatblc 
dctailt bccontc eery apparenu The gen- 
eral underlying thape, however< raaMiai 


eatlly recognlubic, 

aanealatiane arc hercly eiiiblc. 

FIGURE C-1. FRACTIONAL BROWNIAN TRAILE 
(OIMINtlONIO- 'l.llll.Qw l,«au| 



FIGURE C-2. AUTO CORRELATION FUNCTION OF 
generalized brownian MOTION FOR H - 3/4 
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bandwidth. Shannon's source coding theorem states that if the information 
density of the code is greater than the Infonnatlon density of the source, 
the source message can be sent without error. This is an existence theorem. 

In practice, all codes exhibit distortion which rises to high levels as the 
code density approaches the source density. This distortion level can be 
computed from rate distortion theory. 

Now images can set sent with tolerable distortion at a code density, 
corresponding to an Information density of one bit per pixel. Therefore, 
we will assume this image information density in calculating our estimate 
of H. 

Lot the signal spectrum be of the form postulated by Mandelbrot. Let the 
noise be white. Then the signal and noise power spectrums are 


^C(f) - 

lN(r) - N 


(C-1) 


Information rate as a function of the signal and noise spectrums is 
given hy (In* 




df blnlts/second, log.,e ■ 1.443 


(C-2) 


Integrate over G(f) and N(f) to obtain total powers.' 



Pratt, et. al., "Slant Transform Image Coding," pp. 191-209, Fig. 12, 
Digital Image Processing (H. Andrews, Ed.), IEEE Press, 197S. 

Goldman, Information Theory , Sec. 5.2, £q. (43), Prentice-Hall, 1953. 
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Let|M be the ratio of hlgheat to lowest Image harmonic (number of pixels 
per scan 'line). That Is M ■ ^ M > 1. 

SOLVING FOR S AND N 


N 


2iU‘ 


S - 


N 


^2 " ^1 


2U1>, 


2UP. 


-2H --2H 

*^1 ” 


W- 


N 




.-2H 




’N 


(1 - M”^) f. 


N 


<C-A) 


Therefore f 


( 


2HPj, f 


-2H-1 


G(f) 

N(f) 


^ N ^^2 


(2 gRR h\ 

" \ m2H / “ 


ax 


-2H-1 


(C-5) 


where x Is a normalized frequence and " SNR. 

Now for the cases we are Interested In, G(f)/N(f) Is a very large number 
compared to one, except at the band edges near .£2 and. f..' Therefore (C-2) 
can be approximated by the relation 


r ■ e 




■ ef, 


/ 


„ . -2U-1. . 

In (ax ) dx 


An a(x 2 - x^) - (2H + 1) / In x dx 


I 


ef. 


(C-6) 
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The ‘,£2 factor appears because d£ ■ £2 (d£/£ 2 ) " £2 dx. Also when f ■ £21 
X - 1 and when £ ■ £ 1 , x ■ (£,/£„) ■ Carrying out the intogratlon in 

(C-6). 2 



* 

1 

■ '"*2 

(1 - m“^) Ana - (2H + 1) (x An x - x) 

w-1 



M 


(C-7) 

» 


Note that M ^ 1. Also . 


X in X “ 

X-K) 





r - efjCina + 2H + 1) - 2ef (Ana + 2H + 1) (C-8) 

^ » P 


where ■ 2£ . 

2 p 

Now £2 is in cycles/scc (HZ). In order to meet the Nyquist .sampling 
criterion it takes two pixels (a line pair) to de£ine £ 2 * That is £2 ■ 2fp, 
where £p is pixels scanned per second. This observation justifies the second 
right-hand expression in (C-8) . 


Next we must convert information rate to information density, R. The 
relation is 


r(blnit.s/?iec) ^ .. /binits \ 

£ (pixels/scic) " \ pixel J 


So dividing both sides of (C-8) by £^ and using the value of a from (C-5) 


R - 2e 




|An(2SNR) + AnH - 2H An M + 2H + 1 

(C-9) 


] 
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f k 

Negl«ctlng the XnU term, which is small, and solving for H we gat 

(C-10) 

In the case of a typical visual image a representative set of parameters 
could be 

1000 (density range of 3) 

M ■■ 1000 lines (pixels) 

R " 1 bit/pixel 

[h) ■ 

Then 



- 0.*J/|7 -I- y.u(2U00^ 
2t)in(1000) - 1] 


The value of U for the image will be greater than for the object scene 
because the imaging system decreases the high frequency spatial frequency 
content. This Increases redundancy. 

Compensation of U for the Imager modulation transfer function (MTF) will 
be discussed in the next section where H is calculated from measured data. 


C.3 Ul^TEUMXMTION OF U FROM MEASURED DATA 

The only spectrum available were several measured by Charles Hawthorne 
of China Lake. These spcctrums were measured near Faiiioso in the Central 
Valley of California. The terrain type is agricultural. Four fixed beams 
produced resolvable elements which were moved over the terrain by the forward 
motion of the plane. The radiometer signals were digitized and the power 
spcctrums were computed. Figures C-3 and C-'A show results for two separate 
runs. These spectrums were fitted to a power law curves using an BP-97 pro- 
gram. The calculated confidence level was 0.7* This is not too good a number. 
However the dispersion in the data is great. In addition thebe is evidence 
of saturation of the spectrometer channels at low frequencies. 

At any rate. Table C-1 bslsw gives the results of the curve fitting calcu- 
lations as well as the average parameters. The power law curve is of the form 
y ^ ax^-. 

' S'** ■ 
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FIGURE C-3. POWER SPECTRUM VS FREQUENCY 







n>WEfl STECTRUM (CONSTANT X WATTS) 



FIGURE C-4. POWER SPECTRUM VS FREQUENCY 
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. TABLE C-1. POWER DENSITY SPECTRUM FIT TO POWER UW CURVE 


Run 

1 

2 

3 

4 

5 

6 

Average 

a 

24.96 

22.2 

23.62 

19.41 

22.20 

18.06 

21.7 

. b 

-1.86 

-1.86 

-3.05 

-1,42 

-2.18 

-2.79 

2.2 

CuuflXciuat of 
Determination 

0.53 

0.51 

0.64 

O.Bl 

0.82 

0.79 

. 0.7 


The H-value for the image is then 
“lM*GE ■ ^ 


This Is even lower than the value of H computed from the image information 
rate, namely 0.70. 

Uuc to the aparsenesu of the data it is felt that it is reasonable to 
continue to use Mandelbrot value of H, namely 3/4, For the same of compact- 
ness, we will set H + 1 ■ h. Then C(F) ■ F“", 


C.4 SPECTRUM NORMALIZATION 

la ucder tu uhu C(F) tu rupreucat the scene imaged by the radiometer, it 
must be converted to normalized form. In Appendix B it is shown that the 
system responds linearly to radiation intensity. Therefore, the BP power- 
density Spectrum becomes a video voltage spectrum proportional to the RP 
spectrum. The scene spectrum we have been analyzing is given in terms of 
this video power spectrum. Let be this normalized video power spec- 

trum. That is 

G(P)G*(F) « V„(F) V*(F) - S(F) - (C-13) 

0 ' 


The spectrum is two dimensional. The symbol p was chosen to emphasize 
that (C-13) gives the radial component of the spectrum. 


t 
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'i'hUtXalugrai uvur diiu Kpuucturo muwt: be the mean square variance 
scene radiant temperature^ dT<. The rns fluctuations dT vary from 20 
rural areas to 40*K for urban areas* land-water and ice-water scenes 
high contrast with a possible rms values as high as 50*K* 

Thu uonNtant 0 Ih^ found hy Liitegrating the video power spectrum* 
the result equal to AT^ and then solving for 6. 





(p dp d\j<) 


Are AF 


- 1/2 


Since AF’^ - FOV 


AT 


Are FOV 


1/2 


Solving for d 


0 ■ 


AT' 


An 




and 





of 

K for 
exhibit 


setting 


(C-IA) 


(C-15) 


(C-16) 
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OPTIMIZING ANTENNA ILLUMINATION AND BANDWIDTH FOR ' 
MAXIMUM RADTOMRTRTC INFORMATION TRANSFER 

By Jack Gustincic (Aug 1966) 
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iNTROnPCTIOM 

purpose of this tqchnical report is to provide ontenxm desicn 
criteria for an airborne radioaetric sapping systea. The analysis 
employed in detemining optisun antenna illumination and bandvddth 
is considered somewhat novel in that communication theory concepts 
have boon uued for catablishing on objective mathematical evalua- 
tion of system performance. A model of a radiometric mapping sys- 
tem is formulated and the amount of information (in the communica- 
tions sense) available at the system output is calculated as a func- 
tion of various system parameters. The information content is then 
maximized to predict an optimum entenr^ cor>figuration for a receiver 
having a given temperature resolvirg ability. 

Using the information theory approach, it is demonstrated that a 
definite microwave bandwidth and illumination taper must be chosen 
to provide maximum information at the system output. It is showi 
that conventional design approaches which include severe tapering 
of the antenna illumirjation for sideloee reduction and large front 
end bandvridths for increased temperature resolutien do not lead to 
Bftxinun information content. The effects of taper and band'/idth on 
tho apatiiil froeuencics buing scanned by the system are described 
in this analysis and for a tj-pical system, calculations indicate 
that optimum information is achieved *«*ith illuminations near uniform 
(maximum monochrorrntic gain) and band*.ddths around 4^ - Tf’i depending 
on tho resolving capabilities of the receiver and the variance of the 
temperature profile being mapped. 

The use of information content as a measure of system performance 
bos several definite advantages. The system output, for example, 
can be processed in many ways to yield a number of representations 
of tho profile being mapped. A good information content at the out- 
put insures that most processing or coding schemas will have the most 
amount of input data possible and that little data has been destroyed 
Inadvertently by the receiving system before processing. Also, since 
information is an objective o.uantity, it can be maximized without regard 
to subjective quandaries such us hov: such of the main beam is actually 
usefull, etc. 


Clearly, it is impossible to analyse every radiometric mapping system. 
The particular antenna and receiving system discussed in the follovring 
vas chosen bcca'use it was thought to be physically realizable. In any 
ovent the following analysis will serve as an example of haw many simi- 
lar systems' can be treated to provide maximum infozmation transfer. 

i 
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S«ction I beciaa >dth a baaiQ do3oris>tion <?f how microwav® noise from 
the earth is received by the mappini; antenna. 7!ho oonOopts of eoiaaiv- 
Ity end tempovature are introduced and an esear-ession for power received 
at the antenna Is criven. In Soction IX^ this foreuXa is applied to a 
particular nappitur ^stes having on antenna with a separable aperture 
distribution. The antenna pattern is then written as the product of 
pattoma lonsitudinal and tronaverso to the direction of flight of the 
nappinfi: vehicle. 

It ia then shotm, in Eootion III that the t;/o dimensional mappinfi? prob- 
lem reducoo to a sorieo of one dimensional problems involvinjj the trans- 
verse or side-lopkins properties of the antenna alone. In keeping with 
. the state of the art, the trar^svers® ^«tenna conficuration is nssumad 
to bo on array, electronically scanned with ferrite phase shifters. 

Such an ontonna choice allows good pattern control and a fast data gather- 
ing capability. Also, in this section, the concept of spatial frequencies 
ia introduced. Iho temperature profile is eiepandod in terms of its spatial 
fJi*oquonoy coi'iponenta and the spatial frequency tronafor charaotoi'istics 
of the antenna are discussed for the mcnooliromatio case. Finally, the 
effect of microwave bandwidth is calculated and it ia deoons'^ated that 
Incroaoed bandwidth de/7i*ades the fintonna by tapering the apatiai fre- 
quency si»ctrura and introducing unwanted noise. 

In Section IV the effects of the receiver are discussed. A particuler 
receiver cmploitLng straight through detection of the microwave signal 
is aeleoted for the sake of argument. The noise at the output due to 
inadequate temperature resolution is oalculated. The mean square signal 
to noise ratio at the nth spatial frequency 


is then calculated in terms of both antenna end roociver effects. The 
total infornuition at the receiver output is then written as the sun of 
the information oontainod in each spatiol frequency 

% 

I -I ^ *<!-)!] 

n-1 

and this quantity is used as sn objective measure of system performoncoa 
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I^oiaation is plotted as & furstion of bandwidth for various illuiai- 
nation tapers assuninsf a cosine-STuared-on-a-pedeatal family of illumi- 
nation functions. Bandwidth, and tapers which provide naxLnrua informa- 
tion are ppesontsd for a 51 eleDcnt array. 

Zn the last section, the loss of information due to antenna tapering 
ia discussed. It is shown that unwarranted antenna taper equivalently 
increases the receiver noise temperature or increases the microwave 
bandv/idth required to produce a /pLven amount of information at the 
output. A qa?ntity called ''resolution efficiency" is presented as * ^ 
a function of taper. This factor gives a total system efficient 
and is analogous to aperture efficiency in the ease of a monochroma- 
tic antenna. 



lit 
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BASIC COHSIDERATIOltS ' " . , 

The fundumcnUi.X oquationu Ucscribinc the reception of ground mapping 
radiometric information will nov be developed. For purposes of onolysis, 
It will be assxiaed the earth is flat and positioned in the ic-y plane as 
abovn in Figure 1. A receiving aperture « S , is located a height, h, 
above tl^ eart^as shown. A spherical cooruinate system described by the 
ancles Q (h is also indicated. 


it 



Consider an element of the earth, dS, located at the point (x,y). This 
element radiates microwave noise upward into the half space z > 0. The 
noise is produced by thermal effects within the earth and by reflection 
of noise signals frora other sources. It is ass\uned that dS is taken 
large enough so that the noise emitted from it is statistically independent 
of the noise emitted by any other sisiilar element in the x-y plane. The 
alemcnt of dS'is eonpletely described by two additional characteristics; 
aaissivity and temperature. 




I. BASIC COHSIDERATIOUS (Cont>d) 
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Iho tenpcrature of the point (x,y) can he defined in terns of the density 
of power flovinn up and out of <3S. In particular, if dS radiates a total 
average power then the temperature of (x,y) is defined as: 


T(x,y) 


= A. 


2kB 


dP 

dS 


»K 


( 1 ) 


where 


4 


k * Boltznruin ' s constant 
.B ss bandwidth under consideration, Hz 
, A© " center frequency wevelenEth 

Thus , temperature defined in i^Jhis manner' is a measure of the time average 
z component of the Poynting vector averaged over the area dS. 

The emissivity of the element describes any directional radiation charac- 
teristics it might have. The emissivity of dS can be defined as simply 
the antenrm gain Punctioti of the aperture vtiich dC represents. In general 
the emissivity function will be different for the different frequencies 
contained in the spectrum of noise emitted by dS. In particular the follow- 
ing definition for the emissivity of the point (x,y) can be made: 


6(0,^, x,y,CJ) 


power radiated in the direction Q at frequencies 
between LJ end ivJ d LJ " 

power radiated in the direction 0,^'e.t frequencies 
between LJ acid LJ dL/if dS radiated isotropically 
into hnlf space z *>0. 


( 2 ) 


Recognizing that Equation (2) is simply the gain of the aperture dS, the 
noise power received by the mitonnu apertui*e cun cosily be written. Each 
spectral component of noise which leaves dS may be treated as a monochro- 
natic signal transmitted with a gain given by Equation (2). The power, dR j, 
carried away from dS by spectral componeneta of the noise having frequen- ^ 
cios between LJ and LJ * dt*^ enn be calculated assuming the noise follows 
the Rayleigh— Jeans law: 


.. » cu® iU 

where C is a constant. The constant C can be determined by requiring 
that the total power contained in all frequencies be dP, the total time 
average power ^emitted by dB. Thus , 


(3) 


LJ +77b 

P 


dP 


dP. 


u 


fu^ iU Ug stTTio 


W 


U^-TTb 

- 2 - 
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Solvinc Equation (4) tor the constant C the spectral power leaving dS eon 
be written In terinu or the tuin|mruluro as niven by Equation (1). 

dP, . • JJL. T (x,y) iUiS (5) 

^ sTTb TTA" ■ • 

I 


Mow If the power given by Equation (5) leaves dS with a gain given by Equa- 
tion (2) • the power received by the antenna, S. , in the frequency range LJ , 
LJ * dCjis simply: 


dR . 

. — T', A. e ( 0.0.*.v.(j) o( e .<p,u) (6) 

67T 

« I 

* 

where 0( Q %(^%U) is the power gain of the antenna. Substituting for dP 
its value given by Equation (5) and integrating over all frequencies, 
the total time average power received by the antenna fron the.elenent dS 
is then: 


-T(x. y_) 


oo 



U) o(0,p,U)dL/ds 


(T) 


Finally, since dS is statistically independent of position, the tire average 
power received by tho‘ ontennn from nil possible elements of earth is simply 
the Intcgrolof Equation (7) over dS and Is given by: 


rce 


+ 0'0 oo 

( TT ^ fj J 0( g,^,x,y,CJ) C(0,0,Ly) • 


dU/dxdy 


— oo 


( 8 ) 


in which R^ * + y^ and cos Q * h/R, cos 0 ■ x/R sinO, 

Equation (8) describes the auumer in which the antenna receives noise power 
fron the earth and is fundamental to the analysis of a ground mapping 
radiometer. • 


AK A IR BORNE MAPPI MO SY STEM 

Since the width of any airborne antenna is limited by the width of the under- 
side of the aircraft, it seems reasonable to design a radiometric napping 
system with the idea Of enhancing the aperture limited side-looking cnnability. 
Thus, in order to provide pattern control and rapid scon ability, it will be 
assumed that the antenna is chosen to be an electronically scanned array in 
the plane transverse to the direction of the mapping vehicle. Furthermore, 
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AW AiRBOnMK MAPPltig GYGTHM (Cont*d ) 

in order to facilitate the initial formulation of the problem* it will be 
aasiuncd that the emisslvity of the gi^ound is independent of look-angle 
and position. 

For the situation under consideration* scanning effects ore best described 
in terms of the coordinate system shown in Figure 2. A general point (x*y) 
on the earth is located in terms of the direction cosine* u* of the distance 
vector and the longitudinal distance* y* from the airplane. The* coordinates 
(ugy)* are related to the coordinates (x*y), as follows: 

y * y 

u - x7 Vh® + x2 + y2 (9) 

The Jacobian of the transforiaatlon (Equation (9)) cun be calculated directly 
as : ' 





t 


II. AK AIHBOHHB MAPPIUG SYSTKM (Cont^d) 


j 


x,y 




r3 
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10 that the elec»nt of area becomes: 


dx dy ^ 


-f 


du dy 


( 10 ) 




At any sinr;lc frequency, it is assumed that antenna patteni is separable 
and can be written as: 

G((J,0,^). * 0^ (u - Ug) (u,y) JL- (11) 

R 


The function is the transverse gain function of the antenna where Ug 
is the scan produced by an electronic phase shift. is the longitu- 
dinal gain of the antenna and is in general a function of both u and y. 

The factor h/R represents the element factor of the array as it is the 
cosine of the angle the distance vector makes xdLth the antenna normal. 

Substituting Equation (10) and (11) into (8) and Integrating over frequency 
gives the following expression for the power received by the antenna: 




00+1 +0O 



b 




dy du dL/ 



(12) 


The emissivity has been taken as unity end T(u,y) is the temperatu:^e at 
the point (x,y) located by the coordinates (u,y), i.e., T(u,y) is the 

teaperature at the point: < ' 

♦ 

(X . 1L_ [h2 + y2] ^ y) 

1 - 


(13) 


Finally, as the mapping vehicle moves along, the temperature profile 
appears to trwslatc belpv it in the y direction. Tnus, when the airplane 
moves a distance yg from the origin, the received power becomes': 

- 5 - 
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*’pce‘"5>y.* 


IzTT)^ 


CO 4-1 j*?® 

“/ I * 

® «1 • 


ilk) 


CrU.y) du dU 


h® + 


The fom of Equation (l4) indicates that the received power as a function 
Of UgjYg is related to the earth temperature by a linear transfonnation. 
Information about the temperature can be obtained from the received power 
by investigating the manner in which the spacial fz«queneies of the 
temperature ore affected by convolution with the antenna pattern* This 
will be done in the followinc section. 
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RgCKPTION OP SPATIAL PRrl^/UEHCIBS 

9 ho fora of Hquetion (14) Indicatoa that at any elvan position of 
tho aircrnft, y^t tho tranovorno gain function scans an effective 
tonparatura 


T(u) 



^00 

J T Cu, y - yg) Ci (u, y) 

-00 


I? + y^ 


dy 


(15) 


which is the actual teaperature integrated longitudinally with ‘the 
longitudinal gain of the antenna.^ Tho objective nov will be to 
dosign a systou which reproduces ?(u) ns accurately as possible so 
that Equation, (15) can be inverted to obtain the boot possible 
representation of the actual temperature. In this way the total 
two dimensional problem involving u and y first requires, the solution 
to a one dinenoional problem involving u alone. 

la a first step it will be assumed that T(u) is independent of fre- 
quency. The consequences of this assumption will be discussed later. 
Now at any given yg the received power is given by Equations (14) and 
(15) as 

^ +1 

^rcc ^^ 0 ) “ 4^ J / ®t 

0 4 

Ovor the range u ^ 1 the integrated tompprature can bo ozpreased 
in torau of a Eouriur oorioo of opatiol froqueaciea 



If the transverse array is mode up of N 4 - 1 holf-vave spaced* elements 
with frequency independent phase ohiftors^ tho gain function is given 
by 


It can bo shown thot at any particular frequency a half-wave spaced array 
receives all the resolution information which is practicably obtainable* 
Closer or farther spacing produces no additional resolving capability. 
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vhore is the center frequency of the microwave hand* Actually, 

tho c&ixi function is only function u - Ug at the center frequency* 
fact that the spatial frequencies of depend on U/ Uq will 
be shown to stye rise to output degradation* Dhe amplitudes, Gg, 
ere determined by tho excitation of the array. For an array with 
progressive excitation, aj^, K » 1 , 2, • . . N 1, the spatial 
frequency power gain amplitudes are given by 


N + 1 "• n 

•K »K + n 

t 

©u ■ — — — — ■ a ■ 0 , 1 , 

N 1 

I 

K m 1 





(20) 


As will be seen, the Cq will have the effect of weighting the spatial 
frequencies of the temperature and will play an important role in 
determining the total system output. The G|x theh represents the 
spatial frequency transfer characteristics of the antenna at the 
center frequency of the microwave bond. 

Figure shows how antenna illumination effects the Cq for a 51 element 
array. Since Gq is the average value of tho power pattern 

®o "* J ^ (21) 

I 

its value is \mity regardless of the illumination. For a uniform 
array, aK « 1 > the amplitudes decay linearly with increasing n, 


n « 0, 1, • • • N ( 22 } 

, * -8- 



D-12 


ORICINAL PAGE S3 
OP POOR QUAUJY 



Figure 3. S[)atial Frequency Transfer Clwiractcristics for Various Illuminations 












f 
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RrlCEPTIOU OP SPAa?IAL FRStelcrSS (Continued) 

» 

•s ehown by curve (a). Since aXX tho arc aoeuaed to have the 
■aae phase, the area tuidor tho Cq curve ie eiapXy the Baxioun gtkiP 
of tho far field pattern 


■ Cn - <*t (0) - Ot ■(») 


Since a uniform ilXurainatiqn produces oaxioun cain, a tapering of 
tho iXXuminatioii viXX eive a curve having less area than (a). 

Curve (b) illustrates the effect of tapering. An illumination of 
tho form 

•K - .1 + (^’'1 

vas choson for this case. This excitation produces a far field pattern 
having a -41 db side lobe level. As indicated, the tapering tends 
to incroaso tho low spatial frequency responso at the expense of the 
high frequencies. 

Inverse tapering has the opposite effect as ehotm by curve (c)« In 
this case the illumination was taken to bo the inversoly tapered 
function 

ajj - _ 1.143 + ain2 ’ 

vhich produces a 5 db side lobe leve] in tho far field. For the 
inversely tapered illumination the high spatial frequency response 
is increased at the expense of the lows. It is apparent that by 
selecting the array illumination properly, any part of the spatial 
frequency spectrum of Gt (u) can be emphasized but only at the cost 
Of lowering other parts of the spectrum. 

Substituting Equation (l9) into (l6) the expression for the received 

povtr becomes 

» 

N 

J-r... <“.) ■ ^ y" 

n s-=i N 
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RECEPTION 07 SPATIAL FREQUBtlCIRS (OoAtinued) 

It the bracketod exponential In Equation (24) did nG|:'> have to be 
Integrated over the band then Equation (24) would becoae 

i J 

^rec ^^*3) • 27? 

n e - N’ 

♦ ^ 

end the received power would siaply contain the spatial frequencies 
of the teaporature weighted by the amplitudes Gq« However, this is 
not the case. Integrating the exponential gives 

Ms , xJtU- u , nTTu 

f e'"* 0 dtj-sTTse'^ (25) 

Uq-TTs 

where A is the percent bsndvldth 

A = 2TCb/ Uq (26) 

The sin X/X factor in Equation (23) represents an error caused by 
the fact that tho element spacing is not exactly a half wavelength 
over the entire band. This error can be described ma thematically 
by expanding the sin X/X term in a Fourier series as follows 
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III. RSC^JPTIOW 0? SPATIAL PRSCUEacISS (Continued) 
vith ^ 

Si (x) - J ^ it 
0 

m * 

The first three coefficients Co» C) and C2 are shown in Figure 4 
for values of bandwidth ' > 

( 

A<t6/N77' 

For a 50 element array this would correspond to bandwidths up to 
about 10 percents Over this range, the reaaining coefficients C3, 

C4» . . . arc nccligably small. Using Equations (24) and (25) in 
(24) the received power becomes 

+ N 

, ^ V ■ .ATT s»7ru, 

F,oc “ ““ On '0 (n -^ ) 

n - - N 

I* 

N 

+ kB ^ S„[o, +T„.,) 

n « - N 

♦ Cj (a ) (*„ ♦ 2 + *„ . a)] (29) 

where all the Cjj; for K23 have been neglected. 

The effects of system bandwidth ere exhibited in two ways in Equation (29). 
First of all, tho spatial froqusney components in the first sum and 
weighted by the factor Co (nA7r/2J. It is seen from Figure 4 that 
this weighting tends to reduce the amplitude of the higher spatial 
frequencies. The larger tho bandwidth the more pronounced tho 
tepering of tho spoctrua of the Tr, Secondly the second sun in 
Equation (29) represents an addition of noise. The nth frequency 
coaponont of this sum is related to Tjx + l , ^n - 1 » 2 Tjj „ 2* 

Since these aaplitudes are completely independent of T^ they may be 
considered as noise. Figure 4 indicates that the weighting factors 
Cl, C2 and thus the noise, increase with increasing bandwidth. Antenna 
bandwidth then corrupts the system output in two ways: tapering down 

the spatial I'roqueacios of the temperature profile and introducing 
unwanted noia^i in tho output. Further degredetion of the received 
signal by the receiving system will now be discussed. 
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Figure 4. The First Three Coefficients in the Expansion of Sin X/X 







RECSPriON 0? TKB I-^ICRO'.fATC SICT.TAL 


A true moasurexaent of avoraso power reqvdres en indefinitely loizg 
period of tirwj over which instantaneous power is avoragod. Since 
j»:acticcl caesurin^ tices are linited, radioaetric receiving; s^n- 
tems a3J have one factor in cossen* The receiver output varies 
vtth tirae producing; an indication of received microwave power 
which fluctuates about the truo averaso level. In most cases, 
when c^in sjintability is nctrlcctod, the variance of the received , 
power, A corresponding variance in measured teo- 

porature, ®At^, can be estimated by tho formula 

. At® - .(KPT^)® ( 30) 

in which 

K B constant depending on receiver type 

(K » 1,8 for direct detection of microwave Mcnal) 

T^ s= ambient temperature 

r = receiver noise figm*o 

b B post detection band>?idth 

B B predetection bandvddth 

Receiver gain fluctuations also contribute to a variance in the out~ 
put but these effects can be minimised independently by utilizing 
hi^ scan or switching rates. On the other hand, the fluctuations 
described by ( 50 ) arc inversely proportional to the miertr-rave band- 
width and thus represent an isportant consideration in front end 
aystam design. The exact manner in which the receiver output varia- 
tions effect the radicmstric information under investigation is depen- 
dent on tho system configuration, A typical system involving direct 
detection of the microtravo noise will now be discussed. 

The system to be considered is illustrated in Pigure 5« The antenna 
is scanned in ojoicronisa with a storage or display unit at a rate of 
t scans per second. For the purposes of this analysis it will be 
sisxuasd that tho entire 160 dogreo range, -l><.u < 1, is scanned 
although this restriction is not ncccssaxy, ^ 
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IV. REChTTIOIT 0? T?.?, MTCRO’./AVS SIG!IAL (Contiiiuod) 


Array 



Figure 3 


It is assured that f is high enough so that receiver gain fluctua- 
tions only effect the d-c value of the output. As shovna, the micro- 
wave noise from tho antenna is amplified, detected, amplified again 
imd then passed through a processor before storage or display. The 
data processor has boon inclxdcd since it provides an immediate means 
of real time processing at the scan rate'f^. For example, this stage 
could simply bo a filter which ifould variously v?eight the spatial 
frequencies of the input after the scan process converts them to time 
frequencies. 

When the anterna is senuned., tho average povior indicated at the output 
of the video detector begins to vary with time. In particular, tho 
"signal" part of (29) produces a video output tesperatxire. 
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RECSPTl’oN 0? TTra MICROV/AVS SIGH.iL (Continued) 

•fN 

« ttA 3 ^ ^TT f ■ t 

T (t) e • (31) 

. -N 

In this vroy the spatiel frequencies of the input ore converted to 
time froquonciea v.'hich are multiplco of the scan frequency f.^. 

Since the video umplifior must be capable of pessinc all those 
frequencies F the bandvjldth of the anplifier must be chosen as 


b = N f . 


(32) 


Kov; if a total time T is available for transverse scanning there 

will be 7* f total sceno. If in addition, the receiver output is 

averaged over these scans, the variance in the output is ^ven by 

(30) reduced by the factor Mahiivj use of (26) and (32) gives 

s 


At 2 

^ av 


( K P T.)‘ 

TTr 


JL 

B 


( K P 


N 


A 


to) 


Thus the output fluctuation becomes independent of the scan rate, only 
determined by the total nvabor of microwave cycles sampled; Tf^ A« 

Ve are now in a position to examine the information passing through 
ttie system* If the n th frequency component of the "signal” given 
by (31) is sampled it will be found to have a variance 


(n) = 0, ( 


nTTA >2 


) cr , 


a' 


(34) 


where ^T is the variance of the magnitude of the n th spatial 
frequcncy"of the temperature. If the n th frequency component of 
the noise port of (29) is sampled the variance of the noise due to 
antenna bandvridth effects viill bo 






n - 1 
2 


2 2 




n + 2 


•“n - 1 


‘n - 2 


(35) 


« 2 G 


* cr, 2 c,* [ (4^) * ( -^)] 

n ^ ■* 
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IV. R2CSFJI0N 0? TH3 MICHOVJaTO SIorfAL (Continued) 

♦ 

wHore it hao boon nosurjod that ih& vciianco of the tcapcrature 
cooponents docs not vary radically with spatial frec[uoncy* This 
viudonco con bo cntiraatcd by acouainc it io tho swno for all 
spatial fx’oo.uoucios. Thon, if all tho spatial froqucnciea wliich 
could be passed "by the antenna produced a tceicerature profile 
which deviated J T degrees about the avoragOi the voziacce in 
tho n th spatial frequency would be 


cr 2 

n 



(36) 


Finally, assusiing the receiver noise is distributed uniforoly 
about tho spatial frequency cpoctrun, when sampled at the n th 
frequency, the noise ;dll have a variance 





r, ^ "A 


(37) 


Mov/, if, for tho anho of arguaent, the values of the magnitudes 
of each of the spatial frequency components are distributed 
normally among all possible temperature profiles to be invosti- 
catod and if iiio rocoivor noiso is assumed purely random then 
tho veriances given by (34) (35) and (35)‘may be superimposed 
to predict the variance of the receiver output at a given spatid 
froquoncy. Figure 6 illustrates tl'.e oxpoctod vaiuo of the magni- 
tudo of tho frequency spectrum at the output of tho video ampli- 
fier. Tho signal portion of the output spectrum is made up of 
the original spatial frequencies weighted by the gain' and fre- 
quency effects of the antenna. Tho noiso portion of tho spoctrum 
inoludes the effects of antenna bondvddth and receiver temperaturo 
resolution. 
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17. R?.C^yriON OP THS KICR0V;.V»/-£! SIGN.&Ii (Continued) 

In order to obtain an objcctiva ficure of norit for the systeai the 
elasaicol ideas of inforaation theoxy will be invoked end tiM total 
infoiaation present at the output of the video enplifier will be 
defined as 


( 59 ) 



The infornotion contained in the n » o, d - c, coaponent of the out- 
put has been neglected because of receiver instability considerations. 
For a Given receiver and antenna illuaination the inforaation is a 
function of percent ase nicrowave bandwidth, ^ . Substituiting into 
(?9)» the eipreesion for inforaation as a function of bsndvddth becosos 


H 


<a,T 


In 


n a 1 


1 + 




b^(A) + 


■] 


(40) 


whore 


a„(A) 


B.(A) 


Si (-a 


(nTTA )2 


ttA 


)/(-^ 


Si (7T+ JsSiL ) . Si (TT. ) 


t 


+ Si (2rr+ 


nTTA 


)- Si {2 IT - 


nTTA 


K F T. 


JL. 


( 41 ) 
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IV. RSORFTIort OF THB OT.CP.O'./IVS SIGNAL (Continued) 

which io valid for bandvddths up to about 10 percent. The paraceter 
& describes the toaporaturo rosolvins ability of the system relative 
to the temperature profile being investigated. It increases as the 
aq.usre of the receiver noise teaperaturo and inversely as* the variance 
of the temperature profile being scanned end integration tiae. For a 
profile tlmt might have a deviation of / T ■ 60 degrees and radiometer 
figuroo Vi'hich niglit bo typically 

IL B 1.8 

F a 5 db 

7* w 1 ns 

«= IQ 

N = 50 

one firds from (40) that 

R = .004 

Such a system would bo capable of about a 2 degree temperature resolu- 
tion with a band^ddth of about 10 percent, 

Figure 7 sho;;s information as given by (40) plotted for various values 
of tiie resolution factor R, in the ease of a uniformly illuminated array 
^ of 51 elements. The curves show definite morcima v;here the compremiso 
between conflicting antenna and receiver bandwidth effects occur. 

In order to investigate the effects of tapering the illumination of the 
array, a ”cosinc-squared-on-a-pedestal'* distribution function of the 
form 

- B + sin® (42) 

was used to doteraino the snatial freouency power gain aanlitudes, G , 
in (40). Curves of information as a fiuiction of bandiddth and reso- 
lution were obtained similar to those of Figure 6 but for various 
values of tho pedestal heic^it H. The values of pedestal height and 
bandv/idth which nnxinised the inforaatlon were obtained from these 
curves and these results are ahcrum in Figures 8 and 9. 
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Figure 7. Information Vs. Banditddth for a 51 Element lAiifom Array 














Figure 8. IlliBsination Taper as a PunctiOT of Resolution 
Factor for Maxirann Information Transfer 
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Figure 9. ^^ic^qwave Bandidddi Vs. Resolution Factor 
for Maximum Information Transfer 
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IV. 07 *PH3 laCROMm (Conticuad) 

The podcotoX hoifihi; for raEodnvua information io indicated iii Slgure 8 
l 3 i tome of the dh taper from the ccntor to the ends of the array. 

When R ia saoli and tho tyatem is capable of cood tenperaturo resolu- 
tion (lovr receiver noise temperature, Ions intesretion time, etc.) 
the optiRwn b 5 .ndvridth io narrow (Pisura 9; end the opiitoua illumination 
is slichtly inversely tapered. In this case the hifih spatial freg.uen-* 
cios possed by the entonna ore sufficiently above the receiver noise 
fio that they cay be onphasized at tho orpense of tho lot-r freqtuenoics 
to produce a significant increase in inf creation. The ftaphasis is 
achiovod by the inverse antenna taper (hich froq.uonoy vfeinhtinc) and 
, narrow bmidwidth (hi^di fro^iUoncy antenna noise reduction). Por larso 
values of H, fintenna effects are mieaped out by poor noioe integration 
and inforaation is aaxinizod by large bandwidth end positive taper .(lovr 
epatial froquenoy omphaoio). 

It is interesting to note that over a >ddQ range of tho resolution, 
paranoter, R, the illumination for optiinum information transfer does 
not differ from unifora by corn tlum a fevr db. Thus, dccplto sido- 
lobe considerations, anxicua inform tion transfer roughly corresponds 
• to aaxicun'gain. (a. discussion of this remit vrill folloir in the next 
nootion # ) In any oaoe, once the v.nriarce in the ter.'.porature profile 

to bo inva;;ti(;)itod has boon ontinuted, tho ourvoa of Pigurc 8 and 9 
can bo used to select tho opticiuc iUiualnation taper'ond bejidvidth 
for the array. 
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V. COr^ARISON V/ITIt CO:'ITF.:-?OaARY 

The roaulta of tho proviouo analyois havo indicated that luaximun 
InfoTEation tranafer, in the classical sense ^ requires a particu- 
liir choice of aicrovmve bandvridth and antenna taper, dependent upon 
the receiver beinp used and the tenperature profile under investi- 
g.;^tion. These dooic;n criteria are to be contrasted vrith design concepts 
havied on goals v^hich ere aoro subjective than information transfer. 

For example, it seems intuitively plausible that a reduction in 
sidolobe Jovol automatically betters the inforaaticr. gathering capa- 
bilities of the antenna as it rnduccs the amount of spurious po’/;er 
received by the system. Such a statement, hwTever, tacitly assinaos 
that on3.y onergy rocoived by the main beam is of any value in diagnosins 
the target profile. Actually, all the energy at the output of the 
antenna could conceivably be used to describe tho scanned profile. 

In target recognition, for e::ample, tho entire spatial frequency 
spectrum at the system output could be cross-correlated against 
a aynthotic target spectnua using all the available information 
to deduce tho proqon,ts;of the target, not merely that information 
coming from the main beam. , 

When only the oain beam at the antenna pattcro is considered valuable, 
the natural tendency is to choose tlie front end bandwidth as largo 
as possible to provide maxim tamporature resolution in the roceivor. 
Iiarge bandwidths, hovreyer, tend to fill in the sidolobe variations'^ in 
the pattern and destroy some of tho information wliich the antenna 
could obtain v;ith a narrov; bandvddth. 

The actual loss in information due to tapering is shovm in Figure 10 
where the illunination given by (42) was used in (40) to calculate 
mexiraua itifoimiution for vtiriouo pedestal heights. As previously 
oentiojf moxirauci information occurs near a xiniform illunination 
where 1/H =* 0. -As sho:m, tapering in either direction which both 
increases and decreases the sidelobe level, causes a reduction in 
tho maxicnim information which can be passecil by the ^sten. 

Although the numerical value of the information lost nay be open to 
secie subjective interpretation, the objective statement of the infor- 
Bation at the systca output C£m bs used to compare one system against 
another to predict levels of performance. This idea is illustrated 
in Figure 11. V/hcre systems tronaf erring the same amount of informa- 
tion but with different illumination tapers arc Comnriz*od. Thin figure 
dratruitically dononstratoo tlie coot of sidelobe reduction by illuoina- 
' tion tapering. A resolution factor 3. = '*004 vms dioson vihich 

requires a pedestal hoi{^it of H b 4, ,25, -14 db sidelobe level) 

to preduco a maximum information of 55»7 r4cpits at a bandid.dth of 5,5^, 
Now as H is increased end tho illunination becomes more tapered, the 
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V. COI-PAj^ISON WIT? CO:Fi\-!!!?OR.LPY DSSIQIT (Continuod) 

resolution faction of ths syatoin ciuitt be decreased to keep the maxi- 
mum information at 5? #7 nopita. The ratio 

is shovm on curve (b) ao a function of inverse pedestal height and 
is a measure of the resolut,1on efficiency of the syston# Thus, at 
a pedestal height of H « .1, (i/H k 10, -41 db sidolobe level) the 
resolution efficiency is .65 vrtiioh nmans tho intoipration time must 
bo .increased by a factot* of l/.oy « 1.54 pr the rocoiver noiso ficure 
decreased by .94 ^b to maintain fir. output of 55.7 nepits. ?or the 
solce of comparison, curve (a) gives the njonochromatic aj^rtui’e effi- 
ciency and servos to illustruto the roducti.on in monophromatic gain 
due to taper. 

Finally, curve (c) illustrates the reduction in bandwidth vrhich can 
bo m'do 5f tbft opti;-iun taper in used and tJio output Ircpt at 55 *7 
nepits. .As j.nclicatod, tbo system with the pedestal height of .1 
retiuiros about 5 tinea the bsnd^ridth as the system with a pedestal 
hel^it of 4 to produce the sans information at tho system output. 

Tho curvoi) of Pig'.u’c 11 serve to uhovf that sevoro tapering of the 
array ctpiivalently px'oducoa a sicnificant increase in receiver 
noise figure or nnoiutt of nicrewave bcad-.ri.dth needed to produce 
a given amount of output infornation. 
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whio h da p t n dB on t l k » .po a it i «»n »-^ f -t h » tw ^ yi w fa) b ut-not<<)D»*h».4imft4iiflfar»no»-v-. 
It-follows f w itt . ( tOa ^h a t i, l>>l ithin 4^^«>^^»» w lWloy^ ^ «t^B » »o ) e m c n t>^y♦th^oq^ 

..I .— I '1 1 ** 111 '* ' ‘ — 

SO that |/i»| (0-^ |/i»j ><-L ) -i»»pwM> nfa i i U is d s g r— of o e h sw n ee 

eomplex.dtgrtA.of»coJwrs nc a ss>sss ww/ a<><»i^t‘' 4 HadJ ' - x a i s eo l tc d 

^h^ffmUua hin t e m i i yi < 


10.4*3 Calculation of mutual intensity and degree of coherence for light from 
an extended incoherent qua8i<*monochromatic source 
(ft) The Van OiUc.rtr-Ze.mikr. theorem 

Wo Mhull iiuw dolortniiio tlio mutual intensity t/j, and the complex degree of coherence 
fiit for points, and Pg on a screen s/ illuminated by on extended quMi*mono* 
ohromfttio primary source a. For simplioity <r will bo taken to be a portion of a 
pliuiti jMiriillol to.</, mid wo will lUMimiu tliut the intMiiuin iiutwoon the souroo and the 
Hureon is iiomogenuous. Wo also assume that the linear dimensions of a ore small 



Fig. 1U.3. IlUiHtruting tlio van CiTTsuT-ZRitNiiCB thoorom. 


coinpanxi to tlio lilHtaiico 00' Ixitwoon the source and the Horoon (Fig. 10.3), and 
Uiat the angles Ixstwocn 00' and the line joining a typical source point 3 to Pi and Pg 
oro small. 

Imagine the source to be divided into elements dog, dog, . . . centred on points 
<$ 1 , 3p . . ., of linear dimensions small compared to the mean wavelength X. If 
K»i(<) and V^i(t) are the complex disturbances at Pg and Pg due to the element 
d(T«,> ^1^0 disturbances at these points are 

v,(0 » K,(0 - <14) 

Honco , "* 

j{Px,y%) - <y,{i)vm - i<y^i{i)vui)> + m 

m m«i« 
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10.4] 


PARTIALLY OOHSRBNT LIOKT 


000 


Nuw Uu) Uglit vibmliuiiii (iiiiiiiig Iroiii (llltui'uul oluimmlH ui* Uio uuurctt muy bo uuMumud 
to be ite^tioally indepondent (mutually inoohecent}j and of tero mean yalue, >o 
that* 

“ <K^,(0>< W> - 0 when m # n. (10) 

If and Jt^i era the dietancos of P| and Pf from tho aource olemont da^, then 

»'««)- .4. -T-;-. (, --f ) . 

(17) 

where |^m| oharaotorkes the strength and arg tho phase of the radiation from the 
ifith element, t and v is the velocity of light in tho medium between tho source and 
tho eereon. Honoo 


<r..<.)vs.«)> - <^.1 (. - AS (< - %-■)> 


/. m,iW, 


ml 


( 18 ) 


If tho path diftoronno Jt„t — Ii,„i !h Hiuall cornpurtHl to tho cuhorunco lunglli of lliu 
light wo may nogioot tho rotaidution — P„i)jv in tho orgumont of and wo 
obtain from (IS), (IG), and (18) 

•'(^ii^^s) « 2<'^m(0^lS(0> *■ >> ■ — ■ * ( 1 »>) 

Tho quantity <Al„(I)yl;|;(I)> uharautori/.uH tho intonsity of tho radiation from tho 
■ouroo element da^. In. any practical cose the total number of the source elements 
may be oosumed to bo so largo, that wo may regard tho source to Im oiTootively 
miitUiiuous, Donutiiig by i(H) tito iiitoiiHity pur uiiil urua of tho wiuruo, i.u. 
J{Sm)d<Xm — (W) becomesj 

J{P„P,) - J 1[S) -jp- dS, (20) 

li 

whare ilj and J?|denoto tho distances botwoon a typical source point «S and tho points 
Pi and P|i and k ■■ 2wp/v ■■ 2w/A is the wave number in the medium: The complex 
degree of coherence /i(P|,P,) is, according to (20) and (Ob), given by 

I r 

i{Pi ) « ./(P|,/*.) - /(^*i) ’•> •/(/',./',) • 

a c 

are the intensities at Pi and Pg. 


* InMthnmiRft tslwnys iinpUtm n (UKiM^li iMit iiiTctMumrily wi<l») riuiK*% ao«l it(|unti<>ii 

(le) is, in fNel, nul v»li«J fur llMi iiliMiliKiMl <wmi uI Mtrioliy iitUKOiiliniiiintio iigliS. Kor nionooliruiiiatio 
light oiM hsa Kai(i) » LUi< ****^**> ('■sO) ■■ wlioro (/„, sml !/.• era indopondonl ol 

time, so that <P’M(0Vj|^i(l)> <m OLi^ai sad tbia quantity ia in gone^ dilTeraot llrom aeio. 
t Za ganond Am eleo dapenda on diioetion, but for aimpUeity «a naglaet tbia dapondaBoo. 
i Frm now on wa aball ftoquantly uaa tba notation is, dP,. , , for aurfaca alamants oantred 

on tba poiala5, Pi, < 


I a a a ..a 
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510 CRiHoirLBs or orvioa [10.4 

We note that (ho integral ( 21 ) is the Mune m that which oooun in quite a different 
oonneotion; namely in tlie oaloulatlon, on the baaia of the HtnroiH»>!FnK 8 N«< prin< 
oiplo, of tho ooinplox (liHtitrhiMico in tho difTruotion pattern oriaing from diffmotion of n 
Hi^horiiinl wavti on an a|M)i-tnrn hi an 0 |)a<]no Hurcua. Moiu prooiaoly, (21) iinplioa lliul 
iA« compfax degree of cofterenu, whieJt <Uienb*ai the correlation of vibmiion* at a fixed 
point P| and a variable point P| in a plane, illuminated bp an ext e nde d quaei-mono- 
ehrofnatie prinutry eource, is equal to the norntalixed complex amplitude at the eorrte- 
jtondimj point I\ in n cer/ain diffraction pattern, centred on P,, Thie pattern would 
be obtained on replacinq the aouree by a diffracting aperture of the eame eize and ehape ae 
the source, and on filling it with a e^en^ wave converging to P„ the amplitude distri* 
bntinn owr the imm-frimi in the ivjxirlnre /lemg jtroportional to the inteneilg distribution 
ucravi the source, 'L'iila roault svu.h liiiil obtuhiiahi^ by Vak Oittiuix* and kter in a 
aimplor way by Zeomiks.I We ahall refer to it oa the van Oittert-Zemike theorem. 

In moat applications tho intensity I(S) may bo assumed to be independent of the 
position of iS on tho Hurfaco (uniform intonaity), The oorroaponding diffraction 
{iroiduni iH thun lliat of dUtraotiuii of a spliurioal wave of uniform ompUtudo by an 
apurturo of tho samo size and ahape as the aouroe. 

Let (^, 77 ) be the coordinates of a typical source point S, rofeixed to axes at 0, and 
let (X}»T,) and (A'l.Fi) be tho ooordinatos of P| and Pf referred to parallel axes at 
O' (h'ig, hL‘ 1). Thnn, if H ilonoUui Liio iliHtaiivo Of/, 


so that 




li 


I - - 


’ “li/r “ 


m 


Uero only tlio loading terma in Xx/Ji, Yjli, f/P, and ij/P bare been retained. A 
atriotly aimUar oxpreaaion is obtain^ for Ji„ bo thot 


Jti - Jt, 


(A,^-hiV) (X, - y.h 

'2P ~ P 


(23) 


In the denominator of the integrands in (20) and (21), Pi and P, may to a good 
approximation be roplooed by P. Wo also aot 


Tiitm (21) roduiiUN tu 


m 

1 

(Ti - y* 

) 

P 

-p. 


1 

i 

‘4- T|*)-(X,« + 

y**n 

r — 

2P 


e'V 

• 

0 


Ml* 

jjns,vmv 



m 


(24) 

(25) 


(26) 


• . 


Henco if the linear dhnensions of the source and the distanu between Pi and P, ore mall * 
compared to the distance of these points from the sourte, the degree of eokertmee |pu| is 
equal to the absolute value of the normalized Fourier transform of the intensity ^netion 
oftheeouree. 

* P. H. TAK Cirnurr, Phgeiea, 1 (1034), SOI. 
t V. ZKONiail. Physieu. S (1031), 710. 
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Th 0 quantity y> defined by (25) hae a simple interpretation. According to (23) it 
represents the phase difference 2rr(OP| PP|)/^> and may evidently bo nogloctod 
when 

OPj-OP,<I (27) 


For a uniform oircular source of radius p with its centre at 0, (20) gives on iutogra* 
tion (of. § 8.6.2) 


where 


\ V r 


V - ^ I V{Xi - x,f + ( n - r,)* 


V 


2n RXi* -h y,») - (X, 


f 


2H 




(28) 


(29) 


Ji being the BbssiIi function of the first kind and first order.* According to § 8.6.2, 
|2.fi(v)/v| deoreosos steadily from the value unity when v ■■ 0 to the value zero when 
V <m 3*83; thus os the |>oints I\ and P| are aoporated moto and more, the degree of 
coheronco steadily decreases and there is complete incoherence when P^ and P| are 
separated by the distance 

— — .. 0*01 

P^P, » \/(A\ - A,)» -}- ( Ft - Yt)’ = • (30) 


A fu^her inorease in v ro>introducoa a small amount of coherence, but the dogron of 
eohomnoo rttmoiiin smaller tlian <H*I, and i.iiom in rurMier atinpleUi iitoohunmisi for 

V tit 7*02. Siiieo Ji(v) cltangUN sign iiH v puMsuH liiruugh each zen> oiJi(v), the phoso 
Pit "" arg fiit changes there by w; in consequence the position of the bright and dark , 
frmges are interchanged after each disappeoranoe of the fringes. 

The function |2/|(o)/v| dooroasoH steadily from the Vnluo 1 fur v m 0 to 0*88 when 

V ■■ 1, i.e. when 

PiP. = ^~- (31) 


Itogording (t. doxiarture of 12 pur cunt fnnii the kieal value unity as the maximum 
permissible departure, it follows that lAs diameter of the circular area that u illuminated 
almost coherently by a quasi-moruKhromatic, uniform source of angular radi us 
sc ■■ pjR isf 0’16)[/cc. This result is useful in estimating the size of a source needed' in 
experiments on inturfonmoo and dilTraution. 

As an example consider the size of the "area of coheronco'* around an arbitrary 
point on a screen illuminated directly by the sun. The angular diameter 2a which 
the sun’s diso rr>btends on tho surface of the earth is about 0* 32’ 0*018 radian. 

Henoe, if the variation of brightness across tho sun's disc is neglected, tho diameter d 
of the area of cohoronoo is approximately 0‘lGjt/0*000 ISX. Taking tho mean 
wavelength X os 6*6 X 10~* cm this gives d 0*01 mm. 


• No oonftiaion ahould arioo rmin tho r«iut Hint tho nyinhol J iu uIho iiwkI for tho mutuol iuUtnaity, 
so the IsUor always appoon with two sufllxoo or with sovoral argumonta. 

t As early os 1800, E. Vaaen ostiaiatad that the diameter of ths "oirola of ooberaoea" is 
aoBMwfaat omallar than 0*8E)l/s. (Am. SemtiJ. ds rjtfools Normals SupSrisim, t (1180). 291; 
aloe his JCspans d'OpSifiu Phvsiqut (Paria L’Xnprimaria Imperials). 1 (1160), 108.) 
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APPENDIX F 

MAl»PINC GEOMETRY 


System output parameters are .defined in Table F-1, "Mapping Geometry 
Formulas." Inputs parameters and other symbols are defined in the symbol 
list below. 


SYMBOL LIST 








0 

• 

EarCh cancar 

e 

m 

Angla batwaan tila diagonala and 

p 

■ 

RadioBMCar locaclon 



valocity vector 



Location of plxaX arc or edgaa of 
Imiijvr KOV 

N 

m 

Ratio of tlla width and BMixiauB 

<JQ» 

■ 


posalblo width 



Incorsecclon of llnea OP and QQ 

1 

m 

'file ordinal number 

c 

m 



p 

m 

Orbital path langth travaraad par 
laaga buildup 


m 

Tima conatant of phanomanon obaarvad 
by radiomatar 

♦ 

m 

Nadir ongla of boon or of LOS to 
adga of ioagar FOV 

V 

m 

Kudiomotry volocity (input only for ' 
flat-earth caaea) 

a 

m 

Maxiaiia nadir angla 

e 

■a 

Azimuth angle travaraad by acannlng 

beam 

“t 

m 

Earth-concor anglo BubtonilctI by 
ciiu widtli 

All 


AUA<|> • SuUd ungiu uubtundud by 
pixel at radiomatar 

H 


Badlonatar altitude above flat 
earth or orbital alt 

D 

IS 

Scanner antenna width (input only 
for acannara ' 


m 

Earth radlua ■ 6,370 km 

“el 

■ 

Scannur untunnu boigUt (input oniy 

V 

m 

Orbital valocity at aarth aurfaca 


for acannara) 


m 

7.9 Waae 

k 

■ 

Beam width paromatar 

A 

m 

Area of a aap tlla 

K 

ei 

Radiometer performance parameter . 



lifi'fX) ■ First incagar iaea 

than X 
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F.X FU?-KARTH, PENCIL-BEAM SCANNER 

' ‘V.' 

Pigutd F-1 shows Che mopping geomecry in plan and elavaclon for this 
case (Case A). 

From Figure F-1 we see Chat 


FOV - BY 0 


(F-1) 

R - H SEC ^ 


(F-2) 

r - H TAN ♦ 


(F-3) 

-e- 

■ 


(F-4) 

AN - ltA(|. SEC«|< - UA<<» 


(F-5) 

AL - RA<^ - H SEC<^ A<f> 


(F-6) 

'*2 

Aa ■■ Ar AL • H A6 A<f> 

SEC^* - AG SEC^^ 

(F-7) 

t » Ar/V - HA^ SEC^^ V“^ 


(F-8) 



FIGURE F'1, FLAT-EARTH SCANNER MAPPING GEOMETRY 
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M - e/^e 

V ( 

T - t/M - A6 SEC^<^ 6“^ v"^ - — - , |gS. - l 

F.2 PEMCIL-BEAK ORBITAL SCANNER 

The goumetiry for uhXu uuse (Case B) Is presented in Figure F*-2, 
Using triangle OFQ and the lav of sines 

/r + H V , 

\f) » sin”'*’ sin y ■ 6ln”‘*‘(y sin 0 

where ^ has an assigned value and y ■ 1 4- H/r^. 

» 

F sMr/wn 4*^4 mnT* 

'* WM* ««a» 

I 

4> -h (90* 4- n) 4- (90* - I|j) - 180* 

Solving for n 

A 

“1 

^ ■ sin (y sin ^) - 4 

From the #ight triangle GFQ 

« 

R - [H + (1 - COS n) S5C^ ■ 


The ground range to the*plxel arc Is the small circle NQ. Its 




F-5 


(F-9) 

(F-10) 

(F-11) 

(F-12) 

(F-13) 

(F-14) 

length is 
(F-15) 





Then 

s \ 
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1 

Ar ■ 

RA^ SECijj ■ r (1 - cos n) SEC^ SECijJ 

(B 

(F-16) 

AL - 

RA6 ■ ” cos n) A0 SEC^ 

(F-17) 

Aa 

AL Ar - r^(i - cos n)^ A12 SEC^ij) SECt|/ 

(F-18) 

t ■ 

. (UA^. SECiif) *^e s .. 2. 

Ar/V - i ~ “ w (1 - cos n) A*!' SEC^ SEC i|/ 

(F-19) 

M - 

0/AG 

(F-20) 

T ■ 

r ^ 

t/M - ^ (1 - cos n) Afl SEC(J) SEC“*l|( 

(F-21) 


I''. 3 i‘UUiQlJKNCY-i'l-ANU, l-'U^T-iiAimi IMAGliR 


Frequency plane imagers (Case C) measure the two-dimensional spatial fre- 
quency speutrum of the scene. The image is the two-dimensional Fourier trans- 
form of this spectrum. The "snap shot" version measures the spectrum of a 
FOV, forms the image and then repeats. The selected FQV must continue to be 
observed as the Interferometer moves forward one image width. This implies 
the interferomet er an tennas must crack one F5v during this movement and then 
snap to the next FOV. The circular images must overlap to product a swath of 
continuous width x as shown in Figure F-3. 

We wish to produce as great a mapped area per image buildup as possible. 
This mapped area is given by 

i ' 

» 

A ■ (2r sin C)C2r cos C) • 4r^ sin C cos 5 ■ 2r^ sin(24;) (F-22) 


We wish to maximi/.e A, so 

» • 

^ . 2r^ ■" cos (25) - 0 (F-23) 
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In order for the derivative to be zero, ^ must be 45*’. Then sin ^ times 
cos ^ is one h&lf^ So the mapped tile is a square. Its area is 


A - 2r^ CF-24) 


Parameters given for Case A, the flat-earth scanner, in, (F-1) to (F-8) 
apply to Case C as well. However, when these equations are used to describe 
a Fourier imager A6 and A(() are no longer bearowidths. They are orthogonal 
angular resolutions for the final image generated by the Fourier transform 
operation. 

Equations (F-9) and (F'-.IO) arc replaced by expressions for the number of 
spectrum samples which must be taken over the Fourier disk in order to prevent 
aliasing. From Appendix E, Equation (E-26) , this number is 


M - (F-25) 

where D is the diameter of the Fourier disk in wavelengths. 


F-8 
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The >tine per image buildup is 


^ 2H TANi 
^ V 


(F-26) 


So 


£ , 8H TAN4) 

^ ttD^ V 


(F-27) 


also (F-2A) becomes 


A « 2H^ TAN^(|) 


(F-28) 


F.4 FOURIER-TRANSFORM ORBITAL IMAGER (SNAP-SHOT MODE) 
From right triangle OPT in Figure F-4 


*MAX " 




(F-29) 


Let us find the largest angle at earth center which is less than 3,^.y 
and which divides 360® integrally. Half this angle is 


'^’kven 


INT 


UMAX/ 


(F-30) 


Then the angle at 0 subtended by a tile-square width is 


2<j> 


0T - 


EVEN 


N 


(F-31) 


Also, in order to find the fleld-of-vlew we need to know n, half the 
angle at earth-center subtended by the circular image. 


F-9 


FIGURE F-4. FOURIER-TRANSFORM ORBITAL IMAGER 






ORIGINAL PAGE IS 
OF POOR QUALITY 


The itime 

\ I 

available to make 

an image is then 



‘ f 

2r 



r 

e 

N INt/— + ] 
\*MAX/ 

m 


Nr/-iL-\ + i 

V. 


(w) 



(F-35) 


where the value of V comes from Section F.5. 

Since OG ■ r cos n and GQ ■ r sin n triangle GPQ gives 

and the slant range is 


R ■ [H + (1 - cos n) r^3 SEC(|> = Tg(y ~ cos n) 


(F-37) 


The remainder of the output parameters are easily derived. They are 
given In Table F-1. 


F.5 ORBITAL VELOCITY 

For circular orbits » the orbital velocity ratio of two orbits is Inversely 
proportional to the square of the earth radii ratio. Thus: 



(F-38) 


As a special case let one orbit be at the earth surface. Then 




■ V'^ 


(F-39) 


where V ■ 7.9 km/sec and r ■ 6,370 ka^ 
e e 
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F .6 SPINRAD ORBITAL GEOMETRY 

N I 


We see that Art is given by (oiq - where (i) is the satellite orbital 


wu Liiuw xa u/ \ujQ — * 

rate and is the earth angular rate, so 


V t 

V A 

An ■ 2 tt ^ - tat - 2 it — — 

CIR 


(F- 40 ) 


r + H 

CIR ■ C ■ y C - orbit circumference 

r^ e ' e 
c 


(F- 41 ) 


- — ■ 2 ii r^ ■ eurth-secontl pendulum period 
'^e 


tlieu , 


e^ / 2ir \ ^ 

^-v 


n + -jr » 


So from (F^Bd) 


(j) ’ a (^ + A^ 


~1 I 

tan n 
y - cos n' 


(F- 44 ) 



ORIGINAL PAGE IS 
OF POOR QUALITY 


F.7 HP-9,7 MAlPPING GEOMETRY PROGRAM 

Table F-2, which follows^ computes the pertinent geometrical parameters 
for orbital mapping using aperture synthesis or SPINRAD. The effective Inte- 
gration time, teff, for the SPINRAD case Is also found. The first program, 
■RTN/R/S, uses the orbital altitude H to find orbital velocity and period. It 
also finds the nadir angle of the LOS tangent to the earth's surface at the 
horizon, <|imax» corresponding zenith angle of the radius from eart;h 

center to the LOS tangent point, 

The next two programs, LBLA and LBLB, are a pair. Given the nadir angle 
of the edge of the FOV, tf>, they find the earth-center zenith angle or vice 
versa. Tlicy also find the angle of incidence of the llne-of-sight at the edge 
of the FOV, If is known, LBLA enters it into the calculator. * Then when 
LBLB is pressed, n and ij) are found and printed out. Conversely, If n Is known, 
LBLB enters n after it is keyed in and "B" is pressed. Then pressing LBLA 
causes ^ and to be calculated and printed. 

LBLC explores the limits of integration time in the SPINRAD case. A 
sequence of increasing integration times are set up and a set of performance 
parameters are printed for each value of t. First 4 *maX Is printed, the i|)' is 
outputted. This is the nadir angle of original FOV edge after satellite motion 
for timt t. If <j)’ > ♦max* ^ portion of the original terrain disk, whose spatial 
frequency spectrum is being measured, has gone below the horizon. 

The program assumes that the FOV of the linear array normal to its length 
is increased to ensure that radiation from the total terrain disk continues to 
be received after time t. This decreases sensor efficiency, since the wider 
beam receives less power from the disk and -replaces it with power from the 
surrounding area added to accommodate disk motion. This loss is expressed 
as a decrease in Integration time. This decreased time, is computed. 

The satellite motion causes new terrain to replace old terrain at the 
edges of the original disk. These edge effects destroy a certain fraction of 
the image, Ari/r|. This is printed ou^ just before the effective integration 
time, t 

eff 

Finally, the incidence angle of the LOS to the back and front edge of the 
disk will change during measurement. The rearward incidence angle will 
decrease and the forward angle increase. These angle changes, Atli*^ and 
are given. This permits the limitation of t, if these angle changes become 
excessive. 

LBLD is concerned with the "snapshot" mode assumed for pure aperture 
synthesis. The number of square map sectors per orbit is entered. The pro- 
gram then gives the earth-center angle subtended by a map sector (tile) . Next 
it gives the time available to measure the frequency spectrum of the disk 
within which the square is inscribed. Finally it gives the fields-of-view of 
the interferometer antennas. 
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TABLE F-2. APERTURE SYNTHESIS GEOMETRY PROGRAM 
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APPENDIX G 

IMiKKORMANCE Ol<' Al^ElcruiUil-SYNTUESlS IMAGEUS 


G .i IMAGE Si^ECTRUM 

In Section 6, _Part I, of Volume I, "Antennas as Spatial-Frequency Filters," 
it was shown that the KTF of an interferometer with baseline d is 

H(F, d) - Ad"^ ^1 - i. ^ , d-AdiF<d4-Ad 

- 0,0^F£d-Ad OR d + Ad£F<« (G-1) 

When d » Ad (G-1) simplifies to 

F(F. d) - 6(F - d) (G-2) 


In Appendix C, "Microwave Radiometer Image Statistics" it was found that 
the spectrum of the scene expressed as a video power density spectrum is 


S(F) - [— ^1— j p’5/2 _ ep”5/2^ S(F) - G(F) G*(F) (G-3) 

XAtt yFOV/ 


Due to linear response to intensity » the RF power spectrum or video 
voltage spectrum of the radiometer is given by 


Vj(F, d) - H(F, d) G(F) - 6(F - d) G(F) - G(d) (G-4) 


G-2 


ORIGINAL PAGE IS 
OF POOH QUAUTt 


As the baseline length increases £com zero to D, the voltage emerging 
from the radiometer V^(d) Is a perfect replica of the spatial frequency spec-’ 
trum of scene, S(F). Both S(F) and Vx(d) are complex variables. The Fourier 
transform of V(d) is a replica of the projection of the radiant temperature 
distribution of the scene on the baseline. 

The video power density spectrum of the radiometer is then 


1(F) « Vj(F) VJ(F) - H^(F» d) S(F) - 6^(F - D) S(F) - 

(G-5) 


G.2 BIAS ERROR 

In order to find tlie total system error, we need to find the bias error, 
as well as error due to noise. Bias is an error which does not go to zero as 
measurement time becomes infinite. Since the interferometer reproduces G(F) 
perfectly out to F ■ I)* nnd hns zero rcMponso for F > P, the b:las power is 
the LuLegrui oL* tliu power density spectrum from 1) to infinity. 


That is 




(G-6) 


For the sake of compnetness in lutter calculations let 


/rov 



So that 



^ 1 ^ 


- 1/2 



(G-7) 


(G-8) 


^Except for bias due to spectrum curvature, which is ignored for the time 
being. 
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The above calculation asaumet there ie no decorrelation of the algnal 
due to differential tlme-of-fllght (DTOF) of the signal to the baseline ends, 
lit Djut, Huuli duuurirulutXuu uucuca. As the baueline length Increases 6 (p)f 
the spectrum Including the effects of decorrelation p becomes smeller and 
smaller compared to the thre spectrum 6(F) . 


however p If the system bandwidth Is kept small compared to the reciprocal 
of the DTOF, this bias error will be small and can be neglected » Furthermore, 
decorrelation also harms the performance of large flllcd^aperture pencil-beam 
.sennimtK, CaieuJ.atlouH fur tlm effects of lateral decorrelation on a^conven- 
tlonal antenna are given by Beran and Farrent.(^) The magnitudes of the 
effects do not appear to be much different from the eperture-synthesis case. 

Therefore p comparison between the two system types Is not harmed by 
neglecting bias due to decorrelation. 


T.t should be polntod out that the spectrum power reduction vs. spatial 
frequency due to decorrelation Is another MTF. It le multiplied Into other 
transfer functions, such es the MTF due to Image motion, to obtain the system 
M7T. Compensation for some of the performance degradation can be removed dur- 
ing Imnge resLornLlon. Correiutlun bias is treated in more detail in 
Appendix I. 


Tn this nnnivsis, n conm'rvntivi* vnine Tor the syntom 1>;in<lwldth will ho 
chunoii. Thin ImndwidUi will he luudu une-fuurtli the reelprocai DTOF computed 
for the maximum baseline length D and for rays arriving from the edge of the 
FOV, The geometry for finding the DTOF Is shown In Figure G-1 he, low. 



( 1 ) ' 

M. Berran and G. Parrent, Theory of Partial Coherent , Prentice-Hall, 1964, 

Sec. 9.2, pp. 133—137. 
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We see chat 


An n . 

AD D sin 1 1 

-(?) 

for P5V 1 

(G-9) 

D FOV 

DTOF “cl”" 

D FOV 
2f 

f 

» 

• 

(C-10) 


NOTE*. 


C (wavelengths/sec) 


C 

I 


f(HZ) •> operating frequency 


Then 


AL • ' *• 

- 4 DXOF 2D fW 



(G-li) 


U.O NUMUliK OF Sl'ATiAL-FUiiQUliS^CY SAMi'LliS 

This sample number is needed to find the noise variance. For the smooth- 
tnn timo T If) the total ohaorvat ion timo avnllnblis per Iituigu buildup dXvldud 
by the number of samples required. For the time being it is assumed that 
sample density is great enough to meet the Nyquist criterion. Then aliasing 
is prevented . 

( 2 ) 

KO' ' gives the value of this Nyquist sample interval » AF. This value is 


AF - (G-12) 


( 2 ) 3 

' H. KO, "Coherence Theory of Radio-Astronomical Measurements, IE Trans ., 

Vol. AP-15, No. 1, Sec. IV B, pp. 15-16, January 1967. 


page 


2 

The area of a sample cell in the uv-plane ia» therefore, fSV . 

The area of the Fourier disk is it!) 2/4. There the number of spectrum samples 
Tor uniform .snniplu d i ntr Ibut Ion is 



2 2 

nD FOV 

^ 


(G-13) 


This would be the end of the sample number calculation, if we could 
always uniformly sample. However, if the interferometer antennas are moved 
rapidly back and forth along the baseline while the baseline slowly rotates, 
oversampling near the disk center results. This is the situation for radio- 
telescopes which use earth rotation to rotate the baseline. The situation is 
shown in Figure G-2 on the following page. 


The correct number of samples on the half-circle rim in order to meet 
the Nyqulst criterion is (rD/AF) . Each rim sample corresponds to a baseline 
position with (D/dF) samples along its length. So the rotatlng-basellne 
.sample number . Is 


M 


UOT n.'is<0 iiui 


Af 


(n-iA) 


The slowly rotating baseline strategy increases the sample number by four 
to one over the minimum necessary. The best sampling strategy is to rapidly 
rotate the baseline while its length slowly changes. The antennas move with 
constant velocity outward along spiral paths producing uniform sampling. 



FIGURE G-2. GEOMETRY OF SPATIAL-FREQUENCY SPECTRUM 
SAMPLING (SLOWLY ROTATING BASELINE CASE) 
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We are now in a position to write t as 


T 


t ^ 4t 

II FOV^ 


(G-15) 


where t Is the observation time per Image and is number of samples for 
unii'oriii .sampling dunsity. 


C.A NOISF, ERROR 

Next we must trace the generation o£ the video noise power density spec- 
trum. This noise originates in the RF section. Let it have a power spectral 
density of N. This power density is converted to a video power density of N^, 
The total video noise power, alter bandwidth reduction, is dT^. [Remember 
that we suppress (kB) in all expressions! So the true total noise power is 
(k6TB)2.] 

Since there are two antennas and the system is fully coherent K ■■ 1/*^. 
Either the in-phase or quadrature channel (real or Imaginary portion of the 
spectrum sample) will inhibit a noise variance of 


«T^ - ^ 

' 2A1-1 


(G-16) 


The transform process adds the I and Q noise powers to obtain tbo nol.sc 
puwur Lii a pixel. Tliu result is 


5t: 


.J.2 

SYS 

Aft 


(G-17) 


Substituting Af from (C-11) and from (G-15) we get 




"^SYS 


2ft 




(G-18) 
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G.5, DETERMINATION OF OPTIMUM DISK DIAMETER AND MINIMUM TOTAL ERROR 
TuCAl syuttun error iu the eum o£ thermal and bias noise. 

This sum is 


e^(D) - (ST^ + 


(G-19) 


This error must be differentiated, equated to zero and solved for.D in 
order to find the optimum disk diameter. 


d^c^(D) 

dD 


- I + 3f2 


M 0 


(G-20) 


The roHul t l.s 


“ ■ Wi) 


wavuiuugths 


(G-21) 


2 2 

Substituting (G-18) into the expressions for ST^ and b we get 


■ ■ 


n 


l6Tp - *2 


tef’ 


,1/7 ,6/7 ,1/7 
® ^1 ^2 


6-6/7 ,6/7 ,1/7 


CK^) 


- b-'/e CK^) 


(G-22) 


Then the total mean-square error is 


e" - (6^/' + 6-W7) if if . 1 b^ e 




(G-23) 
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APPENDIX H 


l'l•:UI•'OI<MAN(;l•; ok SKUrKS-KKO-AUKAY l*KNCTl.-lJliAM SCANNliKS 


li.i NUiSl:! UlUiOH 

Let the diameter of the circular aperture be D wavelengths. 


A6 


d”^ (rad) 


■ half -power angular resolution 


and 


iLU > AO^ - D ^ (ster) ■■ pencil beam solid angle 
Lut 


*11 

^ i>'UV ' (.sLur) tiulid uii^iu uubLuudud by right 

circular cone scanned by pencil beam 

Tliuii 



JL = IfQV^ 

Ali “ 4 


number of samples per scan cycle 


number of pixels Ln linage 


V - 

H-2 


Then 


(H-1) 


(H-2) 


(U-3) 


And 
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T ■ ^ ■ radiometer integration time (H-4) 


Where t is the image formation time (frame time) . 
Now 


6Tp ■ — (*K^) ■ image noise variance (H-5) 


In the case of aperture synthesis it was noted chat both of the inter- 
ferometer antennas absorbed thermal power within the solid angle H without 
interfering with each other. However, both correlator channels must obtain 
)*owt>r irom tho same pixel in the pencil beam scunner case. Therefore, the 
power per channel in the present case is half the power per channel for aper- 
ture synthesis. The radiometer measures only temperature. That means we are 
using only tJic T-c*l»nnnt‘l output at xcro time delay, l^'or this case K * 1. 

In Appendisc F the question of system bandwidth was discussed. There it 
was pointed out that the bandwidth restrictions due to coherence loss are 
about the same fur both fXlied-aperture arrays and synthetic-aperture cases, 
although some of the details differ. In the fllled-aperture case, decorrela- 
tion uniformly reduces the entire SFD. In the synthetic-aperture case the 
reduction differed for different spatial frequencies. It increased with 
frequency. 

As an illustration of these effects for a filled aperture^ the loss in 
gain vs. coherence length for lateral coherence is shown in Figure H-1. The 
curve comes from IJeran and Par rent’s " Theory oT Partial Coherence " referred 
to in Appendix C. 

We ere interested in longitudinal coherence, not lateral coherence. 
However, it is easy to see how to compute the effect we are concerned with. 

If the antenna is a series-fed array, the signal reaching the output port is 
delayed with respect to the signal entering the array nearest to the port. 
Delay increases with respect to the reference signal as entry into the array 
progresses further and further away from the exit port. 
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RATIO OF APERTURE LENGTH TO COHERENCE LENOTN 

FIGURE H-1. CORRELATION LOSS VS NORMALIZED 
APERTURE SIZE 

I 



Let R(t) be the longitudinal coherence function of the signal, which is 
the transform of the RF spectrum. It may be regarded as part of a signal, 
A^(t), having constant power for all time delays. This signal consists of a 
coherent portion C(t) and a noncoherent portion l(t) . The coherent portion 
Ls of the form R(t). Since power level Is not a function of t, we have 
R(t) + I^(t) ■ A2(t) as shown in Figure U-2. 
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The XncuhurenC signal l(t) will not add amplitude-wise to form an out 
put from the aperture. Its components for differing t*s add power-wise to 
form u noncuUurent result which is scattered out of the aperture and lost. 

The total power in the signal, /RCt) , which adds coherently and is 
received', is 


P 


c 



R(t) dt 


(H-6) 


If there had been no decorrelation , the power received would be 


P 


Reference 



R(0) dt 


(H-7) 


So the loss in gain vs t is 


L - R(t)/R(0) - <t>(t) (H-8a) 

The total fractional loss is therefore 


L 


total 



(p(f) dt 


(H-8b) 


' However, (H-8) is of the general form of the gain loss expressions 
encountered in the aperture synthesis cases. 

So we again will use Equation 11, Appendix VI for system bandwidth. 
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Substituting t in (H-4) and Af in (H-9) into (H-5) with K ■ 1, we get 



(H-10) 


Note that (H-10) is the same as (G-18 ) , Appendix G in the synthetic aper- 
ture case. 


H.2 BIAS ERROH 

In order to determine resolution bias note that D is the cutoff spatial 
frequency of a filled aperture. Then the bias power is 


b 


2 



(p dp d.|;) («K^) 


(H-11) 


This is identical to the integral giving the bias power in the synthetic 
aperture case in (6) , Appendix VZ. 

So ng.nin 


. p,n 


- 1/2 


whore 


AT^ 

Fov 




In addition to (Hffl2) there is bias error because the MTF of a filled aper- 
(iire i.*i iita reuLiuiguiak:. However, some ui* this bias can be removed by 
processing, so it will be neglected. 


H.3 OPTIMUM APliRTUlUi: SIZE MINIMUM TOTAL ElUlOR 

Since the error parametets, f^ and £ 2 * identical with their values 
for the aperture' synthesis case, the optimum D and minimum error expressions 
are the same. also. They. are repeated here for convenience sake. 

I ‘ . 

H 


Kraus, >Radio Astronomy, Sec. 6-9, Eq. 6-65, pg. 170, McGraw-Hill, 1966. 


H-6 
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D 



wavelengths 


’ b^/6 ■ (“K^) 

(6^/' + er^/T) fJ/7 jl/7 (.^2j 


(H-13) 


(^1-14) 


(H-15) 



APPENDIX 1 


RATE>DISTORTION CALCULATION OF ERRORS DUE TO UNDERSAMPLING 
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APPENDIX I 

uA'i:ii-i)isTOia’;i;oN calculation or mu<ous due to undersampling 


i.i mean-square signal error vs information rate 

The mean square error introduced Into a signal as a function of the rate 
of source coding and information transmission through a communications channel 
is best expressed in parametric form. For this mean-square-error criteria and 
for a Gaussian source, this parametric representation is analytically 
tractable.'^’ 


The analysis assumes optimum coding. If coding is optimum, then the dis- 
tortion is an upper bound on the error for aU sources with the same spectrum, 
but with other than Gaussian statistics. 

1.1.1 One-Dimensional Case 


As a warm-up exercise, let us do the one-dimensional case. From Berger, 
Equations (A. 3. 51) and (4.5.52), the parametric representation is 



where a> is radian frequency. 


(I-l) 


(1-2) 


^^^T. Berger, "R(D) for Gaussian Sources," pp. 106-130, Sec. 4.5, Rate 
Distortion Theory , Prentice-Hall, 1971. 

^^^A. Vlterbl and J. Omara, "Continuous Time Sources and Generalizations," 
Sec. 8.4.3, pp. 502-513, Principles of Digital Communication and Coding, 
McGraw-Hill, 1979. 
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liCC the Hpoctruin bo ono»Hldo<l and described by the relation 
♦ (to) ■ w ^ (I"3) 


If there are M lines in the image, we can write 


"max " "“min 


(1-4) 


Furthermore, we Citn uormali^c the frequency range by assuming that 
is unity. This results because this spectrum is self-similar and self-affined. 
(See Appendix C.l.) 

The above considerations permit us to write (Z-l), (1-2) and the Integral 
for the variance power, o^, in the forms. 


4(6) 


■ -/ 


MIN(0, ii»"^) dto - 


2it I 
•'i 


M 


MlN(e, (o“h du» 


/ .I, - M 

MAx(o, in dm - ^ / MAx(o, in 

■F • 


dto 


2 1 r -h . ^ 

L" ■ " 27 / “ 

•'o •'i 


dto 


(1-5) 


The variance, which is used to normalize d, becomes 


y M 

- 1 / -I 


dto ■ - 


to - (h - 1) 
(h - 1) 


M 


\ - .. H ■ („atts) 

27T(h - 1) ^«atts; 


(1-6) 
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c 


Now 6 btcomes greater than wh(^n 


e ■ u 


-h 


, 1 ;i £ H 


(1-7) 


or 


£ 8 5 X 


(1-8) 


So the distortion Integral times 2n becomes 


»-l/h 


h-1 


2 ird(e) 


f” -h -i/h 

I Odd) + / 6) tiu ■ 0(0 - 1) + iji _ 


(1-9) 


Therefore, 


d(e) - 


uu 


h-1 

h 


(h - 1) e ^ M 
27r(h - 1) 


-(U-1) 


(I-IO) 


Normalize the error to obtain 


d(6) 

0* 


h-1 

he ^ - (h - 1 ) 


(I-II) 
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R(0) is evaluated as follows 
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^ni K(0,) 




b)£nu) - talnQ 


, 1 . 



- InO - - h + AnO 

so 


R(0) - ^ nats/UNF 


( 1 - 12 ) 


(1-13) 


where UNF means unit normalized frequency. 


For 6-1 



h - U + 1. - M 


-(h-1) 


1 - M 


-(h-1) 


1 and R(6) 


0 - h h 
Air 


0 


These are the correct limits. Figure I-l compares this case with other 
rate-distortion functions. Quantization also limits information rate as 
show^ in Figure I<^2. 

1.1,2 Two-Dimensional Case 

In the two-dimensional case the differential frequency element in the 
Fourier plane is of the form (pdifidp), where p is radial radian frequency. 


-NATS 



NOTliJ COMI'AHIaON OH MSI* HATE OISTOHTION FUNCTIONS POH 
UNCORHELATEO DATA, NEAREST NEIGHBOR AND MARKOV ■ 

CORRELATION ( r - 0.7B)* AND SPECTRUM 

SCALE APPLIES ONLY TO «1>- SIGNAL SPF.CTRUM CASE 

FIGURE M. ONE-DIMENSIONAL CASES 



FIGURE 1-2. MAXIMUM RATE VS NUMBER OF 
QUANTITiZATION BITS 
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Starting with the Carteal»:n form of the 2-d Integrals from Raft 2, 

Equations (8. A. 45) and (8.4.46), we see that in polar form all Integrals are 
of thu funn. 


I - C 


/ 2ir ^ M r M 

I f(p, e)(p dp d^) - 2irC / f(p, 6) 

A A 


pdp (1-15) 


So we have 


2nd(0) - 0 


/ 


,-1/h 


M 

, ^ , -(h-1) , 

pdp + / p ' dp 


/ • 

A-l/h 


.-1/n 


.-(h-2) 


(h - 2) 


M 


.-1/h 


h-2 

ii o h 
2 ® 


- 0 


h-2 

h 


-JLii 


,h-2 


0 + e - M 


n ,,-(h-2) 


7,11 (h “ 7.) 


(X-. 16 ) 


So 


h-2 


•| 0 ^ + (l - -|) 0 - 


1 - M 


r(h-ii) 


d(0) 


(1-17) 


The variance Is 
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rH 

■ if 


dp - 


2 7j(h - 2) 


M 


1 - 

2w(h - 2) 


So the normalized distortion is 


h-2 


m 

2 

0 


~ 0 + (l - "I) 6 - 


1 - M 


-(h-2) 


M < M < 1 


Next evjilunte tlio niL'u lIiuh 


4iiK(0) 


/ 


,-1/h 


pin ■‘-r— dp. - 


/ 


,-1/h 


(-h plnp - pine) dp 


" It 


(4 - 4 - 


^ >.u0 


,-1/h 


. h 9-2/V+ i 6-*/“ - i 


So 


R(0) 


~ 10s“ ° (nats/UHF) 


I 

1-8 


(1-17) 


<I-18) 


« 

(1-19) 


( 1 - 20 ) 


NINO. OF IMAGE LINESI 
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1.2 HP-97 RATE-DISTORTION FUNCTION PROGRAM DESCRIPTION 

2 

{•'or ulioa(‘ wlUi UP cn leu I a Luc:;, I peugrum Cu uuluuiuLc d(0)/u ' 

•ndi R(6), in the two-dimensional case, below. The symbols used are the same 
as those In use in this appendix. When M Is entered, It Is printed. Then 
*5min ■ is found and printed. K multiplies 0 by a fixed factor, and by 
the operator, on each circuit around the loop. I have been using 2. On each 
circuit 6, d(6)/o^, R(e) and R(e)/R mnv are printed In that order. When 6 > 1 
the progriim stops. The last computation gives R,nax enough to unity to 

permit plotting li(n«x vs uumbor of image linos. The result is shown In Fig- 
ure 1-3. A graph of information rate vs percent distortion is shown in Fig- 
ure 1-4. Ihfonnation rate Is given as a percent of Rnax* ^he parameter char- 
acterizing each curve is number of image lines. 
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LIST OF VICAR IMAGE PROCESSING PROGRAMS 
WITH SHORT PROGRAM DESCRIPTIONS 


From DIGITAL IMAGE PROCESSING, 

Kenneth R. Ceetleman, Prentice-Hall, Inc, 1979 
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IMAGE GENERATION 

Tl)c foljowitiy programi! urc usclut for generating digital images. 

GEN* is a program that produces images in which gray level increases iinesry, 
modulo 256, from the upper ieft>hand comer, at specihed horizontal and vertical 
rntes. 

GRATE generates digital images with verticai or horizontal stripes of specified 
gray levci on a specified background gray level. 

.Sl'O'r* produces tin Image coiilalning u spot. The spot may he Gaussian, con- 
ical, reciprocal, reciprocal squared, exponential, or dome-shaped. Spot location, size, 
and contrast may be specified. 

LOCLI.S* generates images containing 5-by-S*pixcl plus signs at specified loca- 
tions. 

COHER generates sinusoidal coherent noise images and can also add sinusoidal 
patterns to otiicr images. 

. GAUSNOIS'* genurutes random noise Images having a Gaussian histogram of 

specified mean and standard deviation, 

POLYNOIS* generates random noise images having ft histog.ram that is 
SiKciliud hy the user. 

SYNPIC generates an image in complex format having gray level zero. The user 
supplies the real and imaginary values of a few specified pixels. This program can be 
used in conjunction with the inverse i'uuricr transform to produce images containing 
sinusoidal components. 


POINT OPERATIONS 

The following programs implement point operations (gray scale transformations) on 
digital images. 

STRETCH* pcrformi^ general point operations. The transformation may be 
specified os linear, piecewise jinear, or as a cube root or exponential function. It can 
also produce contour lines at specified intervals. 

ASTRTCH2* first computes the gray level histogram of the input image. The 
program then analyzes that histogram to determine the point operation required to 
put the histogram into a specified form. The user may specify either a linear point 
operation or a uniform or Gaussian output histogram. The linear point operation is 
designed to produce a specified amount of saturation at each end of the gray scale. 

MATCH* performs the point operation required to make one image “match” 
another. The point operation may be based upon specified areu common to the two 
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pictures. Altemativelyi for • lingie input imige, the program will make the mean gray 
level of specihed areas conform to specified values. 

ALOIIRAIC OPERATIONS 

The following programs perform pixcl*by>pixcl arithmetic fuiictions. 

PXCAVE* can average up to 10 input images, which may be linearly displaced 
with respect to one another. 

UNITAVIC* divides (he input image into u set of coiiliguous rcciungles and 
averages those rectangles together. In the output image, each rectangle is replaced by 
the average. This program is useful for noise removal from periodic structures. 

DIFFPIC* can add or subtract two images following a linear displacement. A 
specilied linear point operation is also performed on the output image. 

F* is a general'purpose pixel arithmetic routine. The function that relates the 
output image to the two input images has 1 1 specifiable parameters. These may be 
chosen to implement addition, subtraction, multiplication, division, exponentiation, 
and logarithms. The program operates by first generating a 256 x 256 tWo*dimen- 
sionai look-up table, which is subsequently used to produce the output image. 

F2 performs Rcncral arithmetic opera.lioaji on. tme or two iapui inutgeSf The 
arithmetic operation is specilied by a FOKTRAN-iike expression. 

PIXe performs complex arithmetic on two complex-valued digital Images. 

I'iXCRAI) calculates iJic magnitude and angle of the gradient vector of an 
input image. The magnitude is taken as the maximum difference between the current 
pixel and its eight adjacent neighbors. The angle is specified by an integer from 1 
through 8. 

PiXH performs addition, subtraction, multiplication, and division on two input 
images. The output and the two input images may have one or two bytes per pixel. 
The arithmetic is done in a fractional format. 

PIXRMS produces an output image composed of local means or local standard 
deviations of the input image. 

RATIO* is a preprocessor used for comparing two input images. These may be 
compared on the basis of ratio, log ratio, difference, or log differenoe. The program 
generates the parameters required for proper scaling of the output image and fetches 
program F, which performs the specified operation. 

LOCAL OPERATIONS 

* 

The following programs perform a variety of specialized local operations on input 
images. Many of them were designed to remove specific biemishes from images. 
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SAR* can replace speeded reciangutar arcfis with an averaic of surrounding 
gray levels. It is also comntoniy used to copy images from one data set to another. 

AUTOSAK locates pixels that deviate from the average of the pixels above and 
below by niore than a spcciHed tolerance. Such pixels are replaced with the average 
of the pixels above and below. 

QifiAK* can add or siihiract .spccincd values at all pixels within specinedrcctan' 
gulur ureas of the input image. 

PSAR* can add or subtract specified values at alt pixels within specified poly* 
gonai boundaries. 

DKASl^* locates and sets to zero the pixels within connected sets having gray 
level and perimeter below spccillcd thresholds. It is used to remove small objects from 
an image, 

ULDM replaces speciiied ureas in an image by linear interpolation. 

ADL* can add or subtract a specified value to all pixels lying along a straight 
line iKUwcen any I wo pixei.s in the input image. Tiiu diiigonul line is specified by its 
end points. 

ADESPIKE* replaces the gray level of a pixel if it differs from its four nearest 
nqighborK by more than a .specified tolerance. The average of its adjacent neighbors 
is used fur Uie repiaccineni. 


IMAGE MEASUREMENT 

The following programs are concerned with extracting and displaying various mea* 
surements from an input image. 

LIST* cun be used lu list the gray levels or the histogram of an image on a iitke 
printer. 

HISTO generates and plots gray level histograms on a line printer. 

CLSTU* and II1ST2 compute the two>dimcnsionui histogram of a pair of input 
images. 

LPLOT2* produces a graphic plot of the gray levels along a speciiied diagonal 
line in an image. 

LAVE* can average alt the horizontal or vertical lines in an image. The average 
values arc listed on the line printer and output as a onedine image. 

PIXSTAT produces output images which represent the local mean, variance, 
second moment, or standard deviation of an input image. 

LITEXFEK* is designed to calculate the light transfer characteristics using 
input images that arc flat fields at various intensity levels. For specified rectangular 
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•reas within the set of input images, the program computes and plots the mean and 
standard deviation of gray level as a function of the input hrightness Icvei, 

THR£5iHLD* locales and lists areas of an image containing points above a 
specified grr/ ;evei threshold. When such a point is located, a 30 by S0>pixcl area 
surrounding that point is listed on the line printer and set to zero gray level in the 
image. This program is useful for automatically locating small objects, such us stars, 
in an image. 

GRIDLOCA* and GRIDLOCB* are used together to locate the intersections in 
an image containing a rectilinear grid network. These programs arc helpful in the 
geometric calibration of image digitizers. 

DRECK locates line segments in digitized line drawings. 


annotation and dispuy 

Tlte following programs arc concerned with adding various types of annotation to 
images or with circuliiig iniage display. 

ARROW writes arrows into an image at specified locations. 

MARK* superimposes rectangular marks centered at specified coordinates in 
an image. 

SCRIBE* places rectangles around specified areas in an image. 

GRID* overlays a rectilinear grid network into a digital image. 

MAPGRID* overlays an alternating black and white grid (dashed lines) into an 
image. 

OVERLAY* superimposes a latitude-longitude grid onto cartographic projec- 
tion images produced by program MAP2. 

MASK* adds gray scales, pixel coordinate reference marks, label annotation, 
and a gray level histogram to an image in preparation for display. Most of the digital 
images in this book were proiluced using program MA.SK. 

SHADY introduces shading and contour lines into an image. The shaded image 
is actually a partial derivative image taken in a specified direction. This assists visuali- 
zation of slowly varying (low-frequency) images. 

DNSYMBOL replaces each pixel with a square multi-pixel black and white 
* symbol that represents gray level. When the image is displayed, the user Ca4 read the 
gray level ofieach pixel by examining the symbol. 

DISPLAY* produces images by printing on a line printer. Overprinting of 
characters Is used to achieve up to 64 gray levels. 

PRINTPIX also prints images on a line printer. However, it includes facilities 
for demagnification, variable upeet ratio, and a gray scale trantformation. 
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The fotlowing programs are designed for general userspedfled geometric operations. 

CEOM* transforms a specified control grid of contiguous quadrilaterals into a 
specified rectangular grid. This program is particularly efficient when vertical dis< 
placements are small. 

LGEOM* performs the same transformation os GEOM but is more efliicjent 
when vcrlicul dispiaccmcnis arc large. Doth programs use bilinear interpolation of* 
displacement and gray levels in the transformation. 

GEOMA* performs geometric operations where the transformation is specified 
us a mapping of (|uaUrilatcrais into (|uaUriialcruls rather than into rectangles. This 
more general format may alTord a more convenient specification in some cases. For 
example, it is possible to degenerate th$ quudrilatcruis Into triangles and specify the 
geometric transformation ns a mapping between control grids composed of con» 
tiguous triangles, in general. CEOMA runs approximately one*thitd longer than 
GEOM oh comparable transformations, 

POLARECT projects a rectangular Image into a sector of a circle, and vice 

vciiat. 


Rotation and Mannificotion 

Tiic fallowing programs are UcsigneU to perform rotation and magnirication on 
digital images. 

ROTATE* performs 90* clockwise or counterclockwise rotation. 

PLOT* cun perform plus or minus 90* rotation or reflection about the horizontal 
or vertical axis of an image. 

ROTATE2* rotates an image through a specified angle about a specified point 
and places tlie center of rotation at a speuifled point in the output image. This pro* 
gram generates tlie necessary parameters (control grids) and fetches GEOM or 
LCEOM to perform the rotation. 

ANGLE provides (lie .siiecificd rotation and (nmslation necessary to bring two 
images into registration. 

MAG* can generate the parameters necessary to magnify or reduce an image, 
alter its aspect ratio, or skew the top of the image right or left with respect to the 
bottom. Tin’s program fctclics GEOM to perform the operation. 

SIZE* can also magnify or reduce an image or change its aspect ratio. 

RESAMP* reduces the size of an image by r>l;ipping lines add samples. 

APAVG* reduces the size of an image by averaging rectangular arrays of pixels 
to form the new gray level values. 
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EXPAND* enlarges an image by a specified factor N by repeating each pixel in 
an N by array, 

INTERP* magnifies an image by a factor of 2 using bilinear interpolation. 

CLASP changes the linear aspect ratio of an image by magnification in the 
horizontal direction, The program enn use either linear or cubic spline intcrpolutlon. 

Imao* Combination 

The following programs may be used to combine images to form larger ones. 

CONCAT* can combine up to ten input images of the same size. The images 
are concatenated side by side and one above the other. 

INSECT* can combine two images of unequal si'/:c. Specified rectangular areas 
from each of the two input Images arc placed at specified locations in an output 
image of specified size. Where the two areas overlap, the second input image prevails. 
Areas of the output image where neither input image falls are set to 255, 

MOSAIC* is usdul for combining multiple overlapping views of an area into 
one composite output image. The program assumes that the input imases are all the 
same size as the desired output and that the smaller images have already been placed 
in proper position on n hnckground of zero gray level, Where nonzero portions of 
images overlap, the program takes values from the input images in a specified order 
of priority. 

HISTLOC generates an output image from a specified { mensional look- 
up tabic. The gray levels of two input images provide the address otto the table, This 
program is useful for identifying which areas of an ima^ cosrtespond to different 
clusters in a two-dimensional histogram. 

I,YNX joins the top half of one image to the boiitun half of another. 'I'ruiisiution 
it provided to bring overlapping areas into registration. 


Map Projection 

Ti» following programs arc designed for producing map projections of aerial or 
spacecraft images. 

MAP2* is a general cartographic projection program. The program generates 
the necessary parameters to project an image from the camera coordinate system into 
a standard cartographic projection map. The six cartographic projection options are 
Mercator, Lambert conformal conic, oblique orthographic, polar orthographic, 
oblique stereographic, and polar stercographic. The program uses LGEOM to per- 
form the actual projection. 

MAPTRANS can be used to transform a projected image from one cartographic 
projection to another. The program uses LGEOM to perform the necesury geometric 
transformation. 
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MERCATOR generates LGEOM parameters to transform an image into a 
mercator projection, 

PROJECl' generates LGEOM parameters to transform an image from mer? 
cator to orthograpliic projection or vice versa, 

GEOTRAN* performs geometric transformation to effect projection from a 
sphere ontf> a cyiiiiilcr or a plane. It can also perform orthographic projections to 
produce the image that would be obtained if the camera were moved to a specified 
different position. , ' 


CORRELATION AND CONVOLUTION 

The following programs perform tasks that arc implemented with digital correlation 
or uunvulution. 

CROSS* compares specified areas within two input images to determine the 
translation required to register the two images. It computes and prints the sum of 
squared differences between u stationary rectangle in one image and a moving rcc« 
tangle in the second image and the relative displacement that results in mi nim um 
sum>squared difference. The size of the rectangle and the area of search cap be 
specified by the user. It is useful in many ca.scs requiring cross<corrclation orautocor- 
rulutiun, 

REGISTER* generates the GEOM parameters required to bring one picture 
into registration with another. Specified rectangular areas of the images are compared 
by lirsi removing low*frcquency information, then determining the translation that 
results in minimum mean square difference, and finally fitting a polynomial surface 
through the required displacements. A GEOM control grid that spans the image is 
computed from the polynomial surface. 

EILTER* performs general two-dimensional convolution. The input image is 
convolved with a specified rectangular impulse response. This program assumes that 
the impulse response has four-quadrant symmetry. It permits either specified or 
uutumulic scaling of the output gray scale. 

AFILTER^ is similar to FILTER except that no assumptions are made about 
symmetry in the point spread Ametion. 

FILTER2* computes the point spread function of a specified twoHlimensional 
transfer function. The program then fetches FILTER to perform the convolution, 
Thus the user may perform convolution filtering where the filter is spt.clfied in the 
frequency domain. 

GFILTER accepts parameters that describe a twoHlimensional transfer Ametion 
having elliptical cross sections and Gaussian profile. It then fetches FILTER2, which 
computes the corresponding psf and, in turn, fetches FILTER to cany out the con- 
volution. This program makes it convenient to implement lowpass filters of Gaussian 
profile in the frequency domain. 
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UOXFLT* convolves an input image with a f.at*tppped rectangular point 
spread function. The user specifies the size of the p.«;r. Because of the constrained psf, 
this program executes much faster than FILTER. 

BOXFLT2 is similar to BOXFLT except that it can also produce a highpass 
filtered output image. This is obtained by subtracting the lowpass filtered image from 
the input. 

FASTF1L2* can produce a highpass and a lowpass filtered version of an input 
image, The lowpass psf is a two*dimensional rectangular pulse, The highpass filtered 
image is obtained by subtracting a specified fraction of the lowpass filtered image 
from the original. Options to reduce ringing at di.sconiinuilies are included. 

FASTFILI* is a one>dimcnsional (horizontal) version of FASTFIL2. 

« 

MEDIAN produces an output image in which the gray level represents the 
median value of surrounding pixels within the input image. The program is jimited 
to one-dimensional (line-by-line) processing. 

POLYFILT implements convolution with a spatially variant point spread func- 
tion. Up to 10 dUfcrcnl psfs may be supplied. A second input image specifics the 
areas in which each of the psfs are to be applied. 

SMEART^** is designed for Wiener deconvolution of linear motion blur. The 
user specifics the direction and amount of motion blur, and the program calculates 
the psf of the Wiener deconvolution filter assuming white signal and noise. The 
program estimates the signal and noise power spectra from the image and fetches 
AFILTER to implement the convolution. 


FOURIER TRANSFORM COMPUTATION 

Hie following programs perform tasks that involve computation of the Fourier 
transform. 

FFTl* computes the forward and inverse one-dimensional complex Fourier 
transform on a linc-by-linc basis. 

FFITPIX* can be used to display the complex transforms produced by FFTl. 

FFT2*' computes the direct or inverse two-dimensional complex Fourier trans- 
form of a digital image. 

FJ'TIMC* produces a digital image from the complex data set produced by 
FFT2. This permits display of the amplitude and phase of the Fourier transform as 
an image. 

POWER* computes the one-dimensional power spectrum of each line in a 
digital image and displays the rool-mcan-squarc power s^ctrum of all lines on the 
line printer. 

MTF can be used to compute the modulation transfer function of a camera 
system from digitized images of sine wave test targets. 
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FREQUENCY DOMAIN FILTERING 

The following programs implement linear filtering in the frequency domain. 

FREQMOD* multiplies the complex spectrum produced by FFT2 by a user- 
specified transfer function. The resulting complex spectrum can be inverse traiuformed 
by FFT2 to implement linear filtering. The user specifies the profile of an elliptically 
symmetrical, real, nonnegativc transfer function. 

FFTFIT can either multiply a complex spectrum by an input picture or make 
the amplitude of the spectrum proportional to the input picture. Jn either case, phase * 
is unaltered. This program permits frequency domain filtering with unrestricted real 
transfer functions. 

FFTFIL* performs one>dimensional frequency domain filtering on a line-by- 
line basis. Each line is transformed, multiplied by a specified transfer function, and 
inverse transformed. The transfer function may be a bandpass or notch filter with 
user-specified frequency bands. The program can also modify the spectrum by inter- 
polating across specified frequency bands. This program is usefUl for removing 
coherent noise from images. Program POWER may be used to determine the un- 
known noise frequency. 

UNSIIADE* implements highpuss filtering in the frequency domain. It multi- 
plies a complex spectrum (from FFl^) by a transfer function that Is unity everywhere 
except near the vertical and horizontal frequency axes. The transfer function goes to 
zero at the frequency axes with a negative Gaussian profile. This removes low-fire- 
queiicy inrormulioii (shading) from the inverse transformed image. 

SPIKMA5K multiplies the complex spectrum by a transfer function that is 
unity everywhere except in specified small rectangular regions. In thesb regions, it 
takes oil the value zero. This program may be used to remove spikes in the frequency 
domain produced by periodic noise in the image, 

OirrSPlKE'* removes spikes from the spectrum of an image containing per- 
iodic noise. Tliu program first locates spikes in the amplitude spectrum by searching 
for local maxima above u specified severity. It then removes the spike by interpolation 
of the surrounding values. The phase is not altered. 

APODIZE* modifies an image near the edges so as to reduce the effects of 
truncation when the Fourier transform is computed. The program places a quarter 
cycle of a sine function at each end of each line and column of pixels so that making 
the' image periodic in two dimensions does not produce discontinuities. This prevents 
the inlroductioi) of artifact along the frequency axes In the two-dimensional Fourier 
transform, 

OTFl* can compute the horizontal component of a blurring transfer function 
given an image containing a degraded vertical edge. It can also compute the hori- 
zontal component of a transfer function, given the line spread function. The program 
differentiates each line through a degiaded vertical edge to obtain an estimate of the 
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line spread function. These estimates are averaged together and inverse transformed 
to determine a component of the transfer function. 

OTF2* implements two*dimensional frequency domain image restoration by 
Wiener deconvolution. The user supplies the complex spectrum of both the input 
image and the degrading psf. The program assumes white signal and noise spectra, ' 
estimates their amplitude, and computes a Wiener deconvolution transfer function. 

It produces an output spectrum that is the product of the input spectrum and the 
restoration transfer function. The restored image may be obtained by inverse trans* 
forming the output spectrum. 


STEREOMtTRY PROGRAMS 

The following programs arc concerned with stereometry and stereoscopic display, 

RANGE computes a range image, given a stereo pair of input images, 

ELEVMAP generates a topographic map of a .surface, given a stereoscopic pair 
of input images, 

VPROFILE* produces a vertical profile plot of a surface in a range image. 

STEREO* produces a stereo pair, given a brightness image and a range image. 
The user supplies the rigiu eye image and the range image, and the program computes 
the left eye image. 
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I. INTRODUCTION 

This report contains a detailed description and a source listing of 
MEPIX, an iterative Fortran program for maximum entropy Fourier synthesis 
and two-dimensional spectral analysis. The program is distributed to > 
enable Interested researchers to begin exploring two-dimensional maximum 
entropy processing without first undertaking a major software development 
effort. HEPIX was originally designed for use as a radio astronomy 
imaging procedure: it produces a nonnegative image from a finite number 

of noisy, irregularly-spaced samples of an unknown object's Fourier 
transform (the complex visibility function). Such analysis is analogous 
to two-dimensional power spectral estimation from autocorrelation 
measurements. The pfogfam listing in the Appendix Is suff leiently well 
documented that those who wish to use the program for these purposes 
should be able to do so with only cursory reading of the text; however, 
the usual caveats about uninformed use of contributed software are 
operative. While there are no known errors in MEPIX, this does not 
guarantee that there are none, and final responsibility for correct 
operation remains with the user. 

With suitable modification, KEPIX can be used for one-dimensional 
problems of Fourier spectroscopy, spectral analysis of noisy and/or 
irregularly-spaced autocorrelation samples, and estimation of line- 
integrated brightness In radio astronomy. It should also be realized 
that revision of the program to tailor It to special measurement 
geometries may be profitable. If, for example, visibility measurements 
are located at or interpolated onto the vertices of a rectangular lattice, 
the computation strategy can be altered to reduce memory requirements by 
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a factor of about three and execution time by a factor exceeding four. 

The basic reference for this report is •'Maximum Entropy image 
Reconstruct! on by S. J. Wcrnecke and L, R. D'Addario. This paper is 
scheduled to appear in the May 1977 issue of the i^EH Transactions on 
Computers; in the meantime, preprints are available from the authors. 

This paper describes the rationale for and mathematical specification of 
the imaging procedure realized in MEPIX, and it also constitutes one 
chapter of the senior author's dissertation, "Maximum Entropy Techniques 
for Image Reconstruction from interferometer Measurements and Projections 
available from University Microfilms, Ann Arbor, Michigan 1|8106. The 
following section summarizes the mathematical development contained 
therein. 

il. MATHEMATICAL DEVELOPMENT 


Notation 

mj « |m.| exp [j- arg(mj)]; * the Ith complex visibility measurement 

(U|,V|): spatial frequencies of the ith measurement 

M: number of measurements 

2 

0 | ; variance of the ith measurement error, 

a nonnegative, but otherwise arbitrary, image vector whose kth 
component is denoted as fj^ 

f..p; the ideal ‘maximum entropy image. The kth component is 


the image after the n th iteration. The kth component Is 


denoted as f, 


(n) 


N: number of pixels (elements in the image vector) . 

(X|^,y|^): coordinates of the center of the kth pixel 
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Ax^Ay; pixel widths In x** end y-dl rectlonS) respectively 
AA: pixel area (<■ AxAy) 

A: a positive constant selected end perhaps Interactively adjusted 

the user 


Objective Function 

The maximum entropy Image maximizes 

N M 

J(f) - AA E an f - XE -V 

k-1 1-1 0j 

N 


X |mj - AA jJ, f,^ exp[-J 27 r(ujX^ + v,yj]| 


rk " M'k" 


where 


AA E an f, + AA E d, f 


k-1 


k-1 


N N 


k'k 


M 


- ^ E Ep£ - ^E “V im. 


M 


d|^ - 2X E -iyl mjl cos[27r(ujX^ + Vjyj^) arg(mj)3 
I — 1 ^ j 


and 


H 

i?i 27t[u,(x^ - x^) + v,(y^ - 


0) 


( 2 ) 


(3) 


If pixels are arranged on a rectangular grid, there are only about 
2N distinct values In the array and a sum of the form Pjj^ l^f|^ 

for a-1 ,2, . . . Is identified as a two-dimensional convolution sum 
regardless of measurement coverage {(u.,Vj), I - 1,2,. These 
convolution sums are evaluated efficiently by fast Fourier transform 
(FFT) techniques. 


L-6 


ORIQINAI PAOE IS 
OF POOR QUAU7Y 


4 . 


Gradient: 


VJ(f) 


where 


dJ 3J dJ 

.377 577 • • • 177. 


N 


It: ■ 77 * “ S 


ft) 


'i 'I " k-i 
The maximum entropy imege satisfies “ 0. 

Selection of X 

An appropriate choice of X leads to an Image ini satisfying 


M 

I-l cr. 


m 


- M E 


k-1 


4 eI, 


[-J2i(u,*^ ♦ v,y^)] j 


While it is not known how to determine X to satisfy this requirement in 
general, successful use (liDS been made of the formula 


NAA 


2(1 mp/op) 


(5) 


where mp is the measured (or estimated) flux density and Op is Its 
standard deviation. In some cases this value Is too large, and dividing 
it by the number of measurements is helpful . It should be realized that 
iterative adjustment of X may be required to control image-measurement 
compatibility, larger values of X leading to Increased agreement. 

Optimization by Univariate Image Perturbation 

Let be the image after the n th iteration. If only one pixel 
brightness is to be changed in the next iteration, the optimizing step is 
found by setting the corresponding partial derivative equal to zero, the 
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uquoilon c%in bu soivad anolyliGbUy to yield 


r(n+l) 


2Mp 
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( 6 ) 


where Ills the Index of the pixel whose brightness Is changed In the (n-^Ost 
Iteration. For all k i< brightness Is 

changed In each itcrotlon of a univariate .search, It takes at least N 
Iterations to change al rbrlghtnesses. One might set up a correspondence 

f 

between iteration number and adjusted pixel Index such as i ■■ n mod N 4- 1 . 

Equation (6)‘ always yields a positive brightness value since 
p^ ^ Iterative strategy makes J a nondecreasing function of 

Iteration number, Although no rigorous proof Is offered, the univariate 
procedure should converge to the Ideal maximum entropy Image because 
a) the objective function Is unimodal; b) the gradient Is continuous 
(there are no corners In the multidimensional contour diagram of the 
objective function); and c) there Is no artificial damping mechanism 
that might cause Image changes to approach zero prematurely. 

However, there Is no reason to expect this algorithm to be very fast; 
successive steps are perpendicular, which causes a zig-zag trajectory 
Just as with steepest ascent. The univariate method Is also computationally 
undesirable for large N since FFT calculatloit of the convolution sum 
provides no advantage. The reason for this Is that only one element of 
the convolution sum is needed In each Iteration, yet each iteration 
changes all convolution results. 

The attractive aspects of the univariate perturbation scheme are 
automatic satisfaction of the nonnegativity constraint and elimination of 
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tht nted for a search In the per.turbatlon direction (the optimal stepsize 
Is defined inalytlcaUy) , The following multivariate search strategy, 
the one that Is programmed In MEFIX, was developed to retain some of 
the good features of the univariate approach while taking advantage of 
FFT efficiency as well. The flow chart In Fig, I pertains to this 
multivariate method. 


III. THE ITERATIVE ALGORITHM IN MEPIX 


Let be the objective function value and the gradient after 
the n th Iteration. The structure of the multivariate search algorithm In 
MEPIX Is as follows: 


0) Compute and store {d^} (Eq. 2) and {p^ (Eq. 3). 

1) Take the starting Image, , to bo uniformly gray at a level that 
agrees with the measured (or estimated) flux density. Set n <■ 0. 

2) Evaluate and store the convolution sum (using the FFT) 



N 

E 

k-l 


p f^"^ 


for I ■ 1 ,2,, . . ,N. 

3) Compute from Eq. (l). 

A) Calculate VJ^ according tC cq. (A). 

5) Calculate all optimal Univariate moves by a rearrangement of Eq. (6) 


" 23ap 7T " Pi, 


2,2 
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6) Establish a temporary Image a.ccord[ng to 

. ^(n> + ^ |-j(n) . ^(n)] 

"* (n) 

\ 4 heife 0 Is a positive damping factor. If D > 1, all elements of jf' ' 
are guaranteed to be nonnegative. The discussion below contains 
comments on this parameter and on the choice of the fraction R 

introduced In stop 10* 

7) Store the perturbation vector - f and accumulate the 

Inner product 

8. Calculate the convolution sum corresponding to ■. This result can 

overwrite the calculations of step Z. Calculate the objective function 

valued - 

n ■“ 

9. Approximate objective function behavior in the perturbation direction 
by a quadratic fit. To second order, 

J(f^"^ + = J + <VJ ,A^"^ G + i bG^ 

where G is a scalar and b is a curvature parameter whose value is 
estimated by 

b - 2[J^ » - <VJ„,A^'^^> 3 

n n n ™ 

The objective function Is concave so b is always negative in the 
absence of numerical errors* The optimal Image step (relative 

to f^’’^ is specified by 

G* - -<VJ„, A^"^ >/b 

and the corresponding optimal objective function value Is predicted 
to be 
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j* - > ^ 

10) Compute the fractional Improvement (J^ - J^)/(J* - J^). If this 
ratio is greater than the threshold R or if e* ^ 1, make the 
assignments 

f(n+l) . |(n) 



n n + 1 


and go to step if. (In words, make the temporary image the image 
retained after the (n-fl)st iteration.) if the fractional improvement 
is loss than R and c* < 1^ to to step II. 

11) Take the optimal move 

^(n+1) . ^(n) ^ g.# ^(n) 

An alternative computation If step 6 overwrites is 
|:(n+l) _ |(n) ^ ^(n) 

Make the assignment n •4> n •f 1 and go to step 2. 


Discussion 

» 

The current program runs for the number of iterations specified by 
the user, that is, MEPIX does not contain an automatic stopping criterion 
that uses the search history to decide when to terminate. 

The damping factor D (step 6) is incorporated to increase the 
numerical stability of the algorithm. Since the optimal move in the 
perturbation direction Is determined by a simple quadratic 
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(n ) 

approximation, It is Important that the step to the temporary Image jF ' 
be neither too large nor too small relative to changes In the objective 
function value. A fixed value of D can be used, but dynamic adjustment 
of D based on previous experience is expected to give better performance. 
Such an adaptive scheme is available as an option in MEPIX. D Is updated 
according to 


- maxd , (1 - q) min[D^”^//^, (1 + q) 

This formula has evolved from trial and error, but the basic intent Is to 
increase the damping parameter when the move to the temporary Image over-* 
shoots (e* < 1) the optimal Image and to decrease D when the step to 
falls short (e* > 1). The equation ensures that ^ I, which keeps 

Images nonnegative, and it also provides the user with the capability of 
“hard-1 imlting'* the size of the damping parameter change, will 

always be within the range [(1 - q) 0^”^, (I + q) 0^".^ where q is specified 
by the user. If, for example, q ■ .5, the allowed change to 0 Is at most 
50%. If q * 0, 0 remains constant at the original value dictated by the 
user. 

t 

The fractional improvement threshold ratio R (step 10) determines 
when the temporary image is discarded in favor the the optimal move. If 
R ■■ 0, the temporary image is retained whenever it increases the objective 
function value (whenever > J^}. If R ■ I, the temporary image is always 
discarded. R ■ i is recommended— in fact, this value is “hardwired" into 
MEPIX, but it is easily changed — for the following reason: iterations 

taking the optimal step require two evaluations of the convolution sum 
(the primary computational burden of maximum entropy iterations) whereas 
iterations keeping the temporary image require only one and are thus more 
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economic® 1. If the tempo re ry Image almost achieves the optimal objective 
function value In the chosen perturbation direction. It Is wasteful 
spending an extra convolution to ''fine-trim'* the Image when the choice of 
a new perturbation direction will likely lead to an even greater 
Improvement, in practice* this attempt at economy has worked fairly Well} 
savings due to It are on the order of 30t. The key to making this trick 
work Is to select a good ten^porary Image. 

The quadratic approximation of step 9 appears to be excellent If the 
move to the temporary Image Is large enough. When the step size Is tod 
small, the optimal Image Is estimated by extrapolation rather than 
Interpolation, and this leads to greater approximation errors. A sample 
follows: 

n - 23 J - 1.7791 X e*- 0.13765 

D - 5.3 » -1.2796 X 10‘^ J* - 1.8591 x 10“^ 

After taking the optimal move, * 1.8589 x lO which agrees with 
the prediction within two parts In 10^. Even In the early Iterations 
when the trial steps are largest, the approximation Is accurate; 

n - I J - -1.0087 X lo’** c* - 0,kl370 

n 

D - 11.7 J • -1.A562 X lo”^ J* - -5.6508 x lO’^ 

n . 

-5 

After taking the optimal step, - -5.668A x 10 

IV. PROGRAMH I NG DETAILS 

A block diagram of software components In MEPIX is presented In 
Fig. 2. The programming Is modular to facilitate modification. Sub- 
routines In MEPIX are divided Into five functional groups, each of which 

1.-1 A 
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Is described briefly below. Source listings In the Appendix have 
extensive comments so the discussion here Is concerned primarily with 
the Interfaces among software modules and with suggestions for program 
modification. 

* • 
t 

Main Program and I/O Group 

1 

The main program listing In the Appendix indicates the minimum 
structure required for maximum entropy reconstruction. It Is expected 
that new users will want to augment this code for documentation on output 
and error checks bn input or to Integrate the maximum entropy subroutine 
calls within more comprehensive software imaging systems. It might also 
be desirable to permit Image Input for continued Iteration, perhaps with 
a different value of X. , 

There Is little computation within the main program Itself. The 
primary purpose of MAIN Is to reserve storage for arrays, to oversee I/O, 
and to call the subroutines in the overhead group before passing control 
to the main subroutine, ITERAT. Since ITERAT and the subroutines in the 
I/O and overhead groups are called only once, their code could be Inserted 
In-line; however, this has not been done to take advantage of the object- 
time dimension capabi 1 1 ty of Fortran IV. MAIN reserves enough storage to 
accommodate the largest reconstruction task anticipated, but the actual 
dimensions of the Image are part of the Input to the program. This enables 
the user to produce and save a single load module with the flexibility to 
handle different image grid sizes, rectangular as well as square. In a 
multi-user virtual storage environment, this strategy also reduces the 

f 

number of pages- in and pages-out, because array contents are "packed'] and 
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and this in turn reduces the total job time. Two disadvantages of this 
structure are the need for line printer output within subroutine HERAT 
and fairly long argument lists In subroutine calls. 

The two’’dimenslonal Image array Is stored In column- sequential 
format, and its correspondence with the physical x,y-coordlnate system is 
shown In Fig. 3* Note that the outer pixels have their centers shifted 
half a pixel width from the edge of the field of view* The number of 
pixels in the Image and the boundaries of the field of view are Input by 

f 

KAIN. The number of rows and columns need not be equal, but each must be 
a power of two. 

The two subroutines In the I/O group must be supplied by the user. 

i 

PIXOUT outputs the final Image returned by ITERAT. If Intermediate Images 
are needed, they can be obtained by inserting additional calls to PIXOUT 
wilhin ITEKAT. Printout in ITERAT supplies the user wi th Information 

I 

about the numerical behavior of the optimization algorithm. This 

t 

information is described In a later section. 

OATAIN passes visibility measurements and their locations and 
variances to the main program. HEPIX knows that the visibility function 
is Hermitian; therefore, it is not necessary to indicate separate measure- 
ments at conjugate spatial frequencies. MAIN assumes, for purposes of 
computing X(by Eq. 5) and the initial uniformly gray image,- that the first 
visibility corresponds to (u,v) ■ (0,0). This convention is not essential, 
but if it Is to be changed, an additional parameter pointing to the flux 
density measurement is required. Note, a flux density Value (either 
measured or estimated) should always be included to avoid possible 
divergence. No provision is made in the present software for Interactive 
adjustment of X. 
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Overhead Group 

• , 

DXFRM Is responsible for computing the values k ■ l,2,...,N} 
(Eq. 2 ). The implementation of this equation in DXFRM employs some 
trigonometric identities for faster execution. The fact. that the image 
is a rectangular array of pixels enables Use of these identities. No. 
assumptions are made about the locations of visibility measurements.* 
Modification of DXFRM to exploit regularities in measurement coverage 
for increased speed is possible. 

PSF computes the distinct values of '(Pji (Ep* 3/ corresponding to 

- X|^) - n^^Ax 0 1 0 ^ 5. C - l 

(yjj. - V|^) - n^Ay -(R - 1) £ ny < (R - 1) 

where B and C arc, respectively, the numbers of rows and columns in the 
image. Trig identities are again used for speed. The results of these 
calculations are stored in row-sequential format; the first element of 
the output array Is associated with (x^j^ - X|^) ■ 0, (y^ - yj^) ■ (R - 1) Ay 
and the last with (x^ - Xj^) - (c - 1) Ax , (y^^ - yj^) - -(R - l) Ay . 

XFERFN rearranges values computed by PSF and performs a two- 
dimensional FFT of this array. The result is retained for use by 
programs in the high-speed convoiution group. Pre-computing the FFT 
of this array reduces convolution costs by about 3031. 

Two-dimensional FFT Group 

Much of the source code contained In the Appendix Is associated 
with FFT software. The FFT program, CTA (Cooley-Tukey algorithm), is 
transportable only in its Fortran version, which Is significantly slower 

I! ■ 

, ‘f 
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than an assembly language Implementation of the same algorithm. As many 
computer Installations have FFT subroutines that are locally optimized 
for speed, modification of MEPIX to replace CTA with a faster FFT 
realization should be considered. Provided the comments on the CTA source 
listing are studied careful )yi such a substitution ts not difficult. 

It should be noted tha't CTA cannot compute the FFT of a sequence of 
length less than 16, Therefore, use of MEPIX Is restricted to Images 
with at least 8 rows and 8 columns, and the present program cannot be 
applied to one“dlmenslonal problems. Those Involved in one*dlmensIona1 
Fourier analysis are encouraged to rv:ivise the subroutines In the overhead 
and convolution groups. Components In the optimization group need not be 
allorud. 

FFTR arranges for efficient transformation when either Input or output 
data arrays are real. FFT2H Implements the FFT of a two**dimensional real 
array as a sequence of one-dimensional transformations. IFFT2R inverts 
this operation, transforming an array of complex numbers with implicit 
Hermltlan symmetry Into a two-dimensional real array. 


Two-dimensional High-speed Convolution Group 


C0NV2D oversees the high-speed convolution operation. Its basic 
function Is to take the inverse FFT of the product of the output of XFERFN 
and the FFT of the Image. To avoid certain wraparound edge effects, 
padding arrays with zeroes Is required, and the actual transformation 
operations are on arrays with 2R rows and 2C columns. C0NV2D takes care 
of these details automatically. 

MULT20 and CHUL are utility programs. The function CMUL implements 
complex multiplication, and HUt.T2D performs array multi pi Icat ion as required 
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by thg IUgh-sp<s«d convolution tjgorlthm. The storage fornat for data 
processed by FFT2R and IFFT^ZR is rather complicated (see comments 
attached to FFT2R) , which accounts for the length of HULT2D. This data 
structure could be simplified considerably. 

Optimization Group 

« 

The behavior of MEPIX Is fundamentally controlled by ITERAT, the 
subroutIn.o that embodies the Iterative search strategy. Although the 
operations performed by the other subroutines are complex, each operation 
is well defined, and an Improvement that speeds the implementation 
slightly or that cleans Up the code will have relatively little effect on 
the character of the total maximum entropy Imaging paekige. The bulk of 
the source code Is contained In the convolution and overhead groups, yet 

these calculations, assuming they are done correctly, are really 

0 

perfunctory. 

Improvements to the optimization strategy contained In ITERAT can, 
on the other hand, have a profound impact on the program. The iterative 
scheme described In Section III has evolved gradually over the last 18 
months. Experimentation has suggested that this method converges much 
faster than a steepest ascent algorithm, at least’ when applied to radio 

r 

astronomical objects; however, very little is known about why this is so, 
It is expected that substantial improvements to the present algorithm 
await discovery. Thus, the user should not hesitate to Incorporate In 
ITERAT changes that appear desirable. This version of ITERAT has worked 
well for us, but no stronger claims are made. 

OBJFTN and GRAOl are utility routines that evaluate the objective 
function value and gradient, respectively. ROOT calculates the optimal 
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unlvartite moves as described In Step 5 of Section III. JCTREM Is not 
celled In the present version of HERAT. 


Each Iteration results In several lines of printout from HERAT, an 
annotated sample of which Is contained In Fig. A. Much of this Information 
Is of Interest only to those concerned with the behavior of the Iterative 
search algorithm. The most Important parameters for the casual user are 
the fractional Image change, which might be used as a stopping criterion, 
and the rms reconstruction error 



m, 


N 

- AA 2 

k-l 


f|4 expC-J2TT(ujX^ + VjYj^)] 



(7) 


which describes the measurement compatibility of the reconstruction. This 
parameter, which should be approximately one for the final Image, can be 
used to assess the adequacy of. the chosen value of X. if the reconstruction 
accuracy Is not sufficient, the value of X should be increased. Conversely, 
If the reconstruction accuracy Is too good, spurious detail may develop In 
the image, and X should be reduced. 


V. RUN TIME CONSIDERATIONS 

MEPIX has been compiled by the Fortran H-Extended compiler (OPT-2) , 
and a load module is stored on disk. The execution tiine for running this 

« . a 

I 

load module is a function of the number of pixels N/ the number of 
measurements M, and the number' of "'iterations I. The function Is 
approximately 

t - oMN + BNI log2(AN) 

a 

o - l.A X lo“^ sec and 0 - 3.7 x lo“^ sec 

l-22 
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ItERAtlOH? Iteration number 

OBJECT: objective function for current 
Image 

ENTROPY: entropy component of objective 

function 


RHSERR: rms reconstruction error (Eq. 7) 

RMSGRAD; rms gradient 

DAMPING: damping factor 

OBJTEMP; objective function for temporary 
Image 

OBJSTAR: predicted maximum objective func- 

tion In perturbation direction 

FRACSTAR: predicted optimal step size 

GRADPHI: "angle" between gradient and 

perturbation direction 

DXRHS: rms brightness change from current 

Image to new image 


PRMS: rms brightness of new Image 

CHANGE: • ratio of DXRHS to PRHS 

DELPHI: "angle" between current and 

previous perturbation ' 
directions 


KEEP THIS IMAGE means that the temporary 
image Is retained as the result of the 
current Iteration. TAKE OPTIMAL MOVE 
indicates tjuadratic approximation In 
perturbation direction Is used to determine 
new image. 


Fig. i» 
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•re appropriate constants for an IBM 370/168 installation. These 
parameters have been determined experimentally. B is subject to some 
flM/Stuatlon since it depends on the effectiveness of the temporary image 
moves, that Is, how often the temporary Image Is retained. Each optimal 

Step costs an additional FFT convolution. The temporary image was 

* 

retained in 72% of the iterations for the test case that determined the 
cited value of 6. 

With arrays dimensioned to handle ^096 pixels and up to 1000 
visibility measurements, MEPIX needs 36 O K bytes of memory storage. 

This figure could be reduced somewhat by judicious use of equivalence 
statements. In Its most core-efficient realization, MEPlX requires array 
storage to accommodate about IIN floating point values. 
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